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Abstract. We present a global analysis of leptonic and semileptonic kaon decay data, including all recent 
results published by the BNL-E865, KLOE, KTeV, ISTRA+ and NA48 experiments. This analysis, in 
conjunction with precise lattice calculations of the hadronic matrix elements now available, leads to a very 
precise determination of \V US \ and allows us to perform several stringent tests of the Standard Model. 

PACS. 13.20.Eb Decays of K mesons 



1 Introduction 

Within the Standard Model (SM), leptonic and semilep- 
tonic kaon decays can be used to obtain the most accu- 
rate determination of the magnitude of the element V us 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [HE]. 
A detailed analysis of these processes potentially also pro- 
vides stringent constraints on new physics scenarios: while 
within the SM, all d l — > v?lv transitions are ruled by the 
same CKM coupling Vji (satisfying the unitarity condi- 
tion X)fel^jfc| 2 = 1): an d Gf is the same coupling that 
governs muon decay, this is not necessarily true beyond 
the SM. New bounds on violations of CKM unitarity and 
lepton universality and deviations from the V ~ A struc- 
ture translate into significant constraints on various new- 
physics scenarios. Alternately, such tests may eventually 
turn up evidence of new physics. 

In the case of leptonic and semileptonic kaon decays, 
these tests are particularly significant given (i) the large 
amount of data recently collected by several experiments, 
(ii) the substantial progress recently made in evaluating 
the corresponding hadronic matrix elements from lattice 
QCD, and (iii) the precise analytic calculations of ra- 
diative corrections and isospin-breaking effects recently 
performed within chiral perturbation theory (ChPT), the 
low-energy effective theory of QCD. This progress on both 



the experimental and the theoretical sides allows for unique 
tests of the SM that probe very high energy scales. 

An illustration of the importance of semileptonic kaon 
decays in testing the SM is provided by the unitarity re- 
lation 



\Vud\ 2 + \v us \ 2 + \v ub \ 2 



l + A 



CKM- 



(1) 
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Here the Vji are the CKM elements as determined from 
the various d % — > u^ processes, where the value of Gf is 
determined from the muon life time: G M = 1.166371(6) x 
10 _5 GeV -2 [3]- ^ckm parameterizes possible deviations 
from the SM induced by dimension-six operators, con- 
tributing either to muon decay or to d l — > u J transitions. 
As we will show in the following, the present accuracy 
on \V US \ allows us to set bounds on ZIckm around 0.1%, 
which translate into bounds on the effective scale of new 
physics on the order of 10 TeV. 

A detailed analysis of precise tests of the SM with 
leptonic and semileptonic kaon decays has already been 
presented in Ref. 4. However, the significant progress on 
both the experimental and theoretical sides has motivated 
us to perform an updated analysis with three major areas 
of emphasis: (i) the determination of \V US \ from experi- 
mental data, with and without imposing CKM unitarity; 
(ii) the comparison between the values of |V US | obtained 
from data on K — >■ ir(.v (K^) and K — > iiv (K^-i) decays 
and the corresponding constraints on deviations from the 
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V~ A structure of the charged current; (hi) tests of lepton 
universality in K13 decays. 

To carry out this analysis, values are needed for the 
hadronic constants fic/f-K an d /+(0), as discussed in 
Sects. 12.1.11 and 12.2.11 These values are obtained using 
lattice QCD, and various determinations have been per- 
formed. The lattice QCD community, as represented by 
the FlaviaNet Lattice Averaging Group (FLAG) [5HZ], is 
progressing towards convergence on a set of reference val- 
ues, but in the meantime, for the purposes of this work, we 
are led to propose our own. The criteria we have applied 
in averaging lattice QCD results are motivated, but do not 
represent the only set of possible choices. Our adoption of 
these values is intended to illustrate what precision can 
be obtained in testing the SM, given current experimen- 
tal and theoretical results. In particular, wherever possi- 
ble we quote results for quantities such as |T4s|/+(0) or 
IKis/Kidl x fn/fn, which are independent of lattice inputs 
and ready for use as new lattice results become available. 

This paper is organized as follows. The phenomeno- 
logical framework needed to describe Kg 3 and K^ de- 
cays within and beyond the SM is briefly reviewed in 
Sect. [21 The experimental data is reviewed and combined 
in Sect. [3] The results are presented and interpreted in 
Sect. |U 



2 Phenomenological framework 



2.1 Ki2 rates in the Standard Model 



Within the SM, the ratio of photon-inclusive K^ — > t^v 

T £2( 7 )) 



(Ken/ yO to 7r ± — ► i^v (n e2( J decay rates can be written 



as [HUH] 

r, 



\Vus\ 2 P K m K (l~mymlf 



r ne2 ~ \V ud \ 2 f 2 mAl - mj/mir ( + EMJ ' [) 



where /x and f w are the kaon and pion decay con- 
stants, and S EM denotes the effect of long-distance elec- 
tromagnetic corrections. Short-distance radiative effects 
are universal and cancel from the ratio. For pointlike 
kaons and pions, the long-distance electromagnetic cor- 
rections depend only on the particle masses. The dom- 
inant uncertainty on <5 EM comes from terms depending 
on the hadronic structure. Most analyses to date make 
use of the results of Refs. [TO] and [TTJ which were com- 
puted in a model of hadronic structure assuming Breit- 
Wigner form factors for the low-lying vector resonances 
in order to handle the scale matching. These results give 
8 EM = —0.0070(35) (see, e.g., Ref . IT2"|) . Using chiral per- 
turbation theory 13, S] it has been shown that to leading 
nontrivial ordciQ 0(e 2 p 2 ), the structure-dependent correc- 



1 In ChPT, physical amplitudes are systematically expanded 
in powers of the external momenta of pseudo-Goldstone bosons 
(n, K, rj) and quark masses. When including electromagnetic 
corrections, the power counting is in e 2m (p/47r/^-) 2 ". Powers of 
the quark masses count as two powers of the external momenta 
(<D(p 2 ) = 0(m q )). 



tions to <5 EM can be expressed in terms of the electromag- 
netic pion mass splitting. With the relative theoretical un- 
certainty estimated at 25% to account for 0(e 2 p 4 ) effects 
suppressed by chiral power counting, one obtains 



J EM 



-0.0070(18). 



(3) 



With experimental measurements of the inclusive K^ 
and tt£2 decay rates and precise knowledge of the radiative 
corrections, Eq. © can be used to obtain the value of the 
ratio 



V us 



V ud 



Ik 
fl' 



(4) 



2.1.1 Theoretical determination of Jk/J-k 

To experimentally constrain the ratio |V^ s /V^d|, or ulti- 
mately, the value of | V us | itself, a precise estimate of the 
ratio of decay constants Jk/I-r is needed. The analytic 
evaluation of this ratio within ChPT at 0(p 4 ) depends 
on unknown low-energy constants (LECs) and thus can- 
not be predicted with high accuracy. Consequently, pre- 
cise evaluations of fx/f-K are obtained only from lattice 
QCD. However, ChPT still provides useful information 
on this ratio: the SU(3) breaking of Jk/J-k is linear in 
m 2 K — m 2 oc m s — ra u and thus potentially large. For lat- 
tice determinations of fn/fn, this implies that the use of 
very light pions is essential to obtain reliable results. 

During the last few years, new simulations with Np = 
2, N F = 2 + 1, andiV f = 2+1 + 1 flavors of dy- 
namical quarks have been performed by several groups 
using many different lattice QCD formulations, such as 
staggered (MILC Q3]), domain-wall (RBC [IS]), over- 
lap (JLQCD [TO]), and Wilson-like fermions (BMW [H], 
CERN-TOV/CLS Q2], PACS-CS QSJ, and ETMC [201). 
The fundamental characteristic of these new-generation 
unquenched studies is that recent technical and concep- 
tual developments |21| have allowed pion masses well be- 
low 300 MeV to be reached with large physical volumes (L 
up to 4 fm). The PACS-CS [12] collaboration for example 
has already simulated pions as light as m^ = 156 MeV for 
Np = 2 + 1 (degenerate u and d quarks) and clover (first- 
principle lattice QCD) fermions. The resulting PACS-CS 
value, Jk/ fir = 1.189(20) [19 , however, is still plagued 
by large uncertainties due to the small simulated vol- 
ume, Lin^ > 2.3, with corresponding finite-size effects 
<5fse = cxp(— LttIk) w 10%. 

The present status of lattice results for fx/fir [2"2"H3"2l 
[15113 3, 19.34 37] is summarized in Fig.[T] The lightest pion 
mass simulated is listed for each determination, together 
with the smallest lattice spacing used, whenever available. 
The agreement between the different results is remarkable. 
The present overall accuracy is about 1%. 

Among the results for Jk/Itt in Fig. [1] particularly 
noteworthy are the Np = 2 + 1 studies from BMW [37 , 
MILC '09 [S] and HPQCD/UKQCD [3S], for which de- 
tailed evaluation of the systematic errors (chiral, contin- 
uum and infinite-volume limit) have been completed. 
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Fig. 1. Summary of lattice determinations of Jk/ fir 22-32, 15,33, 19,34-37 . The smallest lattice spacing a (when available), 
lightest value of m-n simulated, and type of lattice fermions used are listed in each case All results for which a value of a 
is listed have been extrapolated by the authors to the contiunuum limit. Where two types of lattice fermions are listed (e.g. 
"DWF/Stag"), valence and sea quarks are discretized differently. Unpublished results are indicated by asterisks. 



The recent study of fa/ fir from BMW [37] involves 
simulations of Np = 2 + 1 clover fermions at several vol- 
umes with Lm^ > 4, whereas the continuum limit has 
been estimated by using three different lattice spacings 
(a « 0.065, 0.085, and 0.125 fm). The BMW value [37] is 



value of fa/U [2] is 



fa/U = 1.192(7) stat (6) syst . 



(5) 



An interesting observation from this study is that the 
value of fx/fir a t their lightest pion mass (m^ = 190 MeV 
at a = 0.085 fm) is only 2% below the extrapolated value 
quoted in Eq. ([5]). This strengthens the statement that the 
result is a direct QCD measurement, and not an artificial 
value cured by an ad hoc fit. 

The 2009 MILC result [3T], in contrast to the BMW 
result, was obtained using staggered (AsqTad) fermions. 
Dynamical AsqTad-staggered ensembles have been gener- 
ated by MILC since 2002 [221130] and now represent an 
ample data set for the study of lattice systematics featur- 
ing six different lattice spacings (ranging from a — 0.18 fm 
down to a = 0.045 fm), light pions {m n > 177 MeV), and 
large volumes (Lm^ > 4). This MILC activity marked 
the end of the era of quenched calculations and the start 
of high-precision lattice QCD with light pions. In their 
2009 analysis of fx/fir [2], MILC exploits the new lat- 
tice ensembles with a = 0.06 fm (m^ > 224 MeV) and 
a = 0.045 fm (m, > 324 MeV), together with the older 
a = 0.09 fm (m„ > 180 MeV) data set. The resulting 



fa/U = 1.198(2) stat (±|) syst = 1.197(7), 



(6) 



which is in good agreement with the BMW result (Eq. J5])). 
(We perform the symmetrization of the total uncertain- 
ties in Eq. ([BJ to facilitate error propagation.) However, 
the sizable shift with respect to the earlier 2004 MILC 
result [29 , which made use of a coarser lattice and heav- 
ier pions than were available in 2009, should be noted. 
This is a consequence of the sizable systematic errors af- 
fecting fa //„■: significant shifts in the central values can 
arise when going to lighter pion masses and the contin- 
uum limit (consider also the ETMC results J3HH2Z] in the 
Np = 2 case). For this reason, we consider to be ex- 
ploratory (although important) the studies of fx/ fir in 
which essentially a single lattice spacing is used, includ- 
ing the results shown in Fig. [Q from RBC/UKQCD [33] . 
JLQCD-TWQCD [IB], NPLQCD [35], PACS-CS QI], and 
Ref.[3B] 

Lastly, for the study from HPQCD/UKQCD [35] the 
MILC a = 0.15 fm, 0.12 fm, and 0.09 fm ensembles are 
used; as in the study from MILC, sea quarks are treated 
using AsqTad staggered fermions. However, in contrast to 
the MILC approach, valence quarks are described using 
the Highly-Improved Staggered Quark (HISQ) formula- 
tion. The HISQ QCD action 35 is understood to give 
better lattice scaling — indeed, the MILC Collaboration it- 
self is generating new ensembles with the HISQ action [38] . 



M. Antonelli et al.: Evaluation of \V US \ and Standard Model tests from kaon data 



The f K /f n value from HPQCD/UKQCD [35] is 

SkISt = 1.189(2) stat (7) syst . (7) 

This result is in good agreement with the BMW and MILC 
results. In particular, there is no apparent systematic dif- 
ference between the results obtained using staggered and 
clover fermions. This seems to suggest that possible issues 
associated with the use of staggered fermions (in partic- 
ular, the rooting issue [3HHUJ) are not relevant to the 
determination of fK/fn, at least at the present level of 
accuracy. 

In order to fully exploit the data set in Fig. [TJ we aver- 
age the results of the analyses from BMW, MILC '09, and 
HPQCD/UKQCD discussed above (Eqs. ©, ©, and ©). 
Since these results are consistent, we calculate the average 
weighted by the statistical errors on the individual results. 
This gives our reference central value and its statistical er- 
ror. To obtain the total error on this average, we assume a 
systematic error of 0.006, equal to the smallest systematic 
error quoted among the three inputs. This is justified on 
the basis of the agreement between the results. Adding the 
statistical and systematic errors in quadrature, we obtain 



f K /f„ = 1.193(6), 



(8) 



which is quite consistent with all the results in Fig. [1] in- 
cluding those obtained with staggered fermions and fronr 
preliminary studies. In the above average, possible corre- 
lations between the HPQCD/UKQCD and MILC results 
due to the use of a common ensemble with a — 0.09 fm 
has been neglected. However, since the valence quarks are 
treated differently in these two studies, and the analy- 
ses are completely different, any potential correlations are 
diluted. The above average is consistent with, but has a 
smaller total uncertainty than, both the average from the 
most recent Lattice conference jB], fx/fir — 1-196(10) and 
the preliminary FLAG result [TJ, f K /U = 1.190(10). As 
we note in Sect. [TJ our use of this average to obtain the re- 
sults presented in Sects. ITT51 and B~61 represents a scientific 
choice. Our value for |14 s /V^d| x fii/f-n (Eq. I|55jl) may 
be used with an alternate choice for fx //„. to rederive the 
results of Sects. POI and B~B1 if desired. 

Updates from PACS-CS, RBC/UKQCD, and JLQCD, 
in addition to new results (for example, an Np = 2 + 1 + 1 
result from ETMC [12]), have already been announced 
and will soon improve the present situation. 



2.2 K13 rates in the Standard Model 

The Ku, decays provide ideal channels for the determi- 
nation of \V US \. The starting point of the analysis is the 
expression for the photon-inclusive K — > wiv [Km^ ) de- 
cay rate: 



r K , 3 - 



1927T 3 



k n 2 q 



i + «£S 



(\Vus\ff 

°SV{2)J ) 



(0) I 



A"/' 



(9) 



where Gf is the Fermi constant as determined from muon 
decays, Sew — 1.0232(3) [43ll8] is the short-distance elec- 
troweak correction, Ck is a Clebsch-Gordan coefficient 
(1 for K° and 1/^2 for K± decays), /f°M0) is the 
K° — > ir~ vector form factor at zero momentum trans- 
fer, and Ijct is a phase-space integral that is sensitive to 
the momentum dependence of the form factors. The latter 
describe the hadronic matrix elements 

(Tr(p^)\s^u\K(p K )) = 

{p.+PK)^{t) + (px-p*),*/?"®, (10) 

where t = (pk — Ptt) 2 — {pi + Pv) 2 ■ The vector form fac- 
tor /+(£) represents the P-wave projection of the crossed 
channel matrix element (Q\s^^u\Ktt). The scalar form fac- 
tor fo(t) describes the S-wave projection, and in terms of 
/+(*) and /_(£) reads 



fo(t) = /+(*) + 



t 



-/-(*)• 



(11) 



•k 



By construction, /o(0) = /+(0). Since /+(0) is not di- 
rectly measurable, it is convenient to factor out f*r * (0) 
in Eq. (J3J and then normalize the form factors for all chan- 
nels to f+ * (0), denoted simply as /+(0) in the following. 
The normalized form factors are then defined as 



/+(*) = 



/+(*) 

/+(0)' 



/o(*) = 



fo(t) 
/+(0) 



/+(0) = /o(0) = 1. (12) 



Finally, 



°EM 



long- 
distance EM corrections (Sect. I2.2.2[) and 8^,^ the cor- 



rcprcscnts the channel-dependent 

;Ktt 

SU(2) 



rection for isospin breaking (Sect. l2"T2~3"|) . 

To extract \V US \ from K& decays using Eq. ©, one 
must measure one or more photon-inclusive Kf$ decay 
rates, compute the phase space integrals from form fac- 
tor measurements, and make use of theoretical results for 
/+(0), #em, an£ l ^suT2V ^ e discuss the evaluation of these 
different ingredients in the following. 



2.2.1 Theoretical determination of /+(0) 

The vector form factor at zero momentum transfer /+(0) 
is the most critical hadronic quantity required for the de- 
termination of \V US \ from K(3 decays via Eq. ([9]). By con- 
struction, /+(0) is defined in the absence of electromag- 
netic corrections. More explicitly, /+(0) is defined by the 
K° —> 7r~ matrix element of the vector current, Eq. ([TO"]). 
keeping kaon and pion masses at their physical values. 

In this section, we restrict our discussion to the eval- 
uation of /+(0) in the isospin limhn (m u = md). This 
hadronic quantity cannot be computed in perturbative 
QCD, but is highly constrained by SU(3) and chiral sym- 
metry. In the chiral limit and, more generally, in the SU(3) 



The choice of the K — > it form factor as the common 
normalization is motivated by its smoothness in the m u — md 
limit (see Sect. 12. 2. 3|) . 
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limit (m u = m,d = m s ) the conservation of the vector cur- 
rent implies /+(0) = 1. Expanding around the chiral limit 
in powers of light quark masses one can write 



/+(0) = l + / 2 + / 4 + 



(13) 



where /„ = C(m™ d >s /(47r/,r) n ), and f 2 and / 4 are the 
next-to-leading order (NLO) and next-to-next-to-leading 
order (NNLO) corrections in ChPT. The Ademollo-Gatto 
theorem implies that [/+ (0) — 1] is at least of second order 
in the breaking of SU(3) or in the expansion in powers 
of m s — m, where m = (m u + md)j2. This also implies 
that f-z is free of 0(p A ) countcrterms in ChPT and can be 
computed with high accuracy: f% = —0.0227 [441I45] . 
Difficulties arise with the calculation of the quantity 

Z\/ = / 4 + / 6 + ... = / + (0)-(l + / 2 ), (14) 

which depends on the LECs of ChPT. The original quark- 
model estimate of Leutwyler and Roos [35] gives Af — 
-0.016(8) and /+(0) = 0.961(8). More recently, analytical 
calculations have been performed to evaluate the NNLO 
term / 4 , writing it as 



h = A{n) + f 4 \ loc (n), 



(15) 



where A{fi) is the loop contribution, with /j, the renormal- 
ization scale, computed in Ref. SSJ and / 4 | loc (^) is the 
0(p 6 ) local contribution involving 0(p 4 ) and unknown 
0(p 6 ) LECs. To estimate the latter term, various ap- 
proaches have been used, including a quark model [46] . dis- 
persion relations [37], and 1/Nc estimates [3BJI3H]- These 
studies obtain the results Af = 0.001(10), -0.003(11), 
0.007(12), and 0.009(9), respectively. Af is found to be 
compatible with zero from these studies. Relative to the 
quark-model estimate from Ref. [35] these new results are 
obtained using more sophisticated techniques and feature 
better control of systematic uncertainties. The resulting 
values for /+(0) are summarized in Fig. [2 As can be seen, 
the uncertainties are no smaller that those from the orig- 
inal estimate of Ref. 45 s which illustrates the difficulty 
of calculating Af to below the 1% level using analytical 
methods only. 

Figure[5]also summarizes present results for /+(0) from 
lattice QCD [50H56] . As can be seen in Fig. [5] lattice 
results, while in remarkable agreement with the original 
quark-model estimate of Ref. 1451 give somewhat smaller 
results for /+(0) than do recent analytical calculations. 

On the lattice, /+(0) is determined using three-point 
functions at non- vanishing momenta, whereas fx/fn is de- 
termined from two-point functions at rest. Because of this 
additional difficulty, current lattice calculations of /+(0) 
are less mature than those of /#-/ '/„•. In particular, most 
results shown in Fig. [5] were obtained with only one lattice 
spacing and with heavy pion masses. Additionally, only 
one calculation of /+(0) exists with Np = 2 + 1: that from 
RBC/UKQCD [551155]. While BMW, MILC, and HPQCD 
currently have interesting results for fx/ fn, these groups 
do not yet have results for /+(0). 

Nevertheless, the special chiral properties of /+(0) 
make it possible to obtain lattice estimates with rela- 
tive uncertainties comparable to those for fx/ fir- Among 



the lattice results for /+(0) in Fig. [2] particularly note- 
worthy are the studies from RBC/UKQCD '07 [55], 
RBC/UKQCD 10 [IS], and ETMC [53]. The correspond- 
ing values of /+ (0) are respectively 

/+(0) = 0.9644(33) stat (34) syst _ oxtr (14) syst _i at , (16a) 

/+(0) = 0.9599(33) stat (+^) syst _ cxtr (i4) syst . lat , (16b) 

/+(0) = 0.9560(57) stat (62) syst . (16c) 

(The two contributions to the systematic errors on the 
RBC/UKQCD results are from extrapolation uncertain- 
ties and lattice effects, respectively.) In Ref. [55] and in 
their update [56], RBC/UKQCD make use of a simula- 
tion with Np = 2 + 1, but a rather coarse lattice spacing 
(a = 0.11 fm) and a lightest pion mass of m^ « 330 MeV. 
Even though the use of smaller lattice spacings would be 
advisable, the corresponding error on the SU(3) breaking 
of /+ (0) seems to be under control to within the stated sys- 
tematic uncertainty (see, e.g., discussion in Ref. [5]). A criti- 
cal issue is the chiral extrapolation from their points with 
ra w > 330 MeV. In the new study from RBC/UKQCD 
[56j . there is an attempt to quantify the systematic error 
from the extrapolation, but a better understanding of the 
NNLO terms in the chiral expansion^ as well as additional 
Np = 2 + 1 simulation points at smaller pion masses, are 
needed to fully address this delicate point. 

The recent ETMC study [53] makes use of very light 
pions, as well as different lattice spacings and volumes. 
Both SU(2) and SU(3) chiral fits are investigated and 
give compatible results, which is in contrast to the find- 
ings of RBC/UKQCD, in which the SU(2) chiral extrap- 
olation for /+(0) has few points with m^ < 400 MeV 
and looks problematic [55] . In summary, the ETMC result 
features a more thorough exploration of systematics than 
the RBC/UKQCD result. However, the ETMC simulation 
still has Np = 2, and the final uncertainty on /+(0) is in- 
flated to take into account unknown systematics related 
to the quenching of the strange quark. At present, many 
other groups are occupied with finalizing their studies of 
fa/fir- As attention returns to /+(0), further progress on 
understanding these systematics should be possible. 

For the numerical analysis in Section [3] we use as our 
reference number 



/+(0) = 0.959(5), 



(17) 



which is our symmetrization of the recent RBC/UKQCD 
result |56| . However, we note that this value is fairly rep- 
resentative of the results and spread of values in Eqs. (fT5|) . 
This number is also basically consistent with the average 
of the N F = 2 ETMC and Np = 2 + 1 RBC/UKQCD 
results as quoted at the most recent Lattice conference 
/+(0) = 0.962(5) 0. 



2.2.2 Electromagnetic effects in Kgz decays 

The results of the most recent calculation [58] of the long- 
distance electromagnetic corrections to the fully-inclusive 



3 Ref. [57] describes exploratory NNLO fits using 
RBC/UKQCD results for f (t). 
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Fig. 2. Present determinations of /+(0) = f+ * (0) from analytical or semi-analytical approaches 
|50H56| . Unpublished results are indicated by asterisks. 



and lattice QCD 



Mode 


«£& (%) 


K% 


0.495 ± 0.110 


K*k 


0.050 ± 0.125 


Ks 


0.700 ± 0.110 


K% 


0.008 ± 0.125 



and strong isospin-breaking tests (as in Sects . l4~2l and 14 . 3|) : 



Table 1. Electromagnetic corrections to the fully- inclusive 
i^3( 7 ) rate [58] . 



i-Q 3 ( 7 ) decay rates for each of the four modes (the 5^m 
in Eq. ©) are listed in Table [TJ These values were ob- 
tained to leading nontrivial order in chiral effective theory, 
working with a fully-inclusive prescription of real photon 
emission. For the low-energy electromagnetic couplings 
appearing in the structure-dependent contributions, the 
recent determinations of Refs. [5H] and [BD] were used. The 
errors in Table [T] are estimates of higher-order contribu- 
tions that are only partially known. The associated corre- 
lation matrix was found to be 158' 



-1.000 +0.081 +0.685 -0.147^ 
+ 1.000-0.147+0.764 
+1.000 +0.081 
+1.000, 



(18) 



It is also useful to list the uncertainties on the linear com- 
binations of S^m that are relevant for lepton-universality 



°EM 


°EM 


°EM 


-6 K±f * 

"EM 


°EM 


°EM 


8 K±tl 

"EM 


"EM 



= (-0.205 + 0.085)%, 
= (0.042 + 0.085)%, 
= (-0.445 + 0.160)%, 
= (-0.692 + 0.160)%. 



(19a) 
(19b) 
(19c) 
(19d) 
the 



The corresponding electromagnetic corrections to 
Dalitz plot densities can be found in Ref. [551 It is im- 
portant to note that the corrections to the Dalitz distri- 
butions can be locally large (up to ~ 10%), with con- 
siderable cancellations in the integrated electromagnetic 
corrections. Hence a proper implementation of the electro- 
magnetic corrections in the analysis of the experimental 
data is essential, in particular for a reliable extraction of 
the form factor parameters. 



2.2.3 Isospin-breaking corrections in K(% decays 

In Eq. ©, the same quantity / + (0) = /f O7r ~(0) (form 
factor at zero momentum transfer) is factored out for all 
decay channels. The isospin-breaking corrections are then 

where 



included via the term containing S^J, 2 -. 



*K°*- 

°SU(2) 



= 0, 



°SU(2) 



/r^°(Q) 



- 1. 



(20) 
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This term can be related to the ir°-r] mixing J5TJE5] At 
leading order (0{p 2 )) [44] 



°SU(2) 



|i withi? 



while at NLO in the chiral expansion (0(p )) 



m d - m v 

4\ 



x K+k° 

°SU(2) 



3 1 



--(1 + X/ + 4m+OK)) 



4i? 



(21) 



(22) 



where \p A ~ 0.219 is an 0(p A ) correction calculable in 
ChPT [22]. Am is a correction (starting at 0{m q )) to the 



ratio m K /m^: 



'K 



1 / m< 

. o 1 + — 

m~ 2 \ m 

where Q 2 — (mf - - - 2 
Eq. (1221) can be rewritten 

3 J_ 

4Q2 



(l + A 



M 



91 

R 



(l + A 



M) 



rh 2 )/(m 2 d — m 2 ). Using Eq. 



?K + +> 
°SU(2) 



m 



K 



2 



in 



(23) 

423), 

(24) 



which shows how <5suf2") ' s essentially determined by the 
double ratio Q 2 (the dependence on m s /rh is suppressed 
by the smallness of Xp 4 )- One can extract Q 2 from the 
analysis of the decay r/ — > 3ir or from the kaon mass split- 
ting. A recent analysis using the latter method gives [2H] 

Q = 20.7 ±1.2, (25) 

and thus (using m s /m = 24.7 ± 1.1 and including 
0(e 2 p 2 ,p 2 ) corrections [5T] to Eq. (IMjl ) 



SU(2) 



= 0.029 ±0.004. 



(26) 



This is based on an evaluation of the low-energy electro- 
magnetic couplings [55] leading to a large deviation of 
Dashen's limit [53]. Note that previous analyses of r\ — > 3ir 
decays [64, give higher results for Q, and hence central val- 
ues for 8g V ,2\ below the lower edge of the range of values 
quoted in Eq. (j2"5|) . New analyses of this decay based on 
recent data [65] are in progress [55; 67_ and should shed 
light on this issue. 

As a final note, the precision reached in the measure- 
ment of the Kgs decay rates and in the determination of 
the corrections 5^ allow iSfyA to ^ e determined directly 

By means of Eqs. 

J SU(2) 

be used to derive interesting constraints on the quark mass 
ratios. 



from data, as discussed in Sect. 

and (|2~4"1) . the empirical determination of #£ T A can then 



2.2.4 Parameterization of the form factors 

The last ingredient for the determination of | V^, s | from 
Eq. ^j is the calculation of the phase space integrals, Iki 

(27) 



x [f + it) + 



3™^k - 2l 



(2t + mj)\ 



fv/oW 



with A K tt = m 2 K - ml, A = [t - (m K + m^) 2 ] [t - (tjik — 
m,r) ], and £ max = (itik — rn^) 2 . In order to calculate the 
integrals, knowledge is required of the normalized vector 
and scalar form factors defined in Eq. (TT21 . The form fac- 
tors can be determined from fits to the measured distri- 
butions of the K13 decays in t or some equivalent variable 
using a given parameterization for the form factors. 

Among the different parameterizations proposed in the 
literature, one can distinguish two classes [55] • Parameteri- 
zations based on a systematic mathematical expansion are 
to date the most widely used. In this class (Class II by the 
nomenclature of Ref. 68), one finds the Taylor expansion 






+ -X'" — 



(28) 



where A' + and A" are the slope and the curvature of 
the form factors, respectively. Another Class-II parame- 
terization is the so-called z-parameterization of Ref. [551 

In Class-II parameterizations, the parameters describ- 
ing the higher order terms of the form factor expansion 
are free to be determined from data. In practice, this addi- 
tional freedom greatly complicates the use of such param- 
eterizations. As noted in Ref. [Till if a quadratic parameter- 
ization is used for both the vector and scalar terms, fits to 
experimental data will provide no sensitivity to Aq because 
of the strong parameter correlations, especially between 
Aq and A '. For this reason, existing power-series fits use a 
parameterization in A' + , A" , and Aq (sec Eq. (|28p). It has 
been shown in Ref. 71 that in order to describe the form 
factor shapes accurately in the physical region, one has to 
go at least up to the second order in the Taylor expansion. 
This is quantified in Ref. 1701 if the same K^ spectrum is 
fitted using both the linear (Ao) and quadratic (A , Aq) 
parameterizations, one typically finds Aq ~ Aq + 3A '. Ig- 
noring the quadratic term increases the phase space inte- 
gral by about 0.15%. In addition, as discussed below and 
in Sect. 13.5.31 for tests of low-energy dynamics involving 
the Callan-Treiman theorem, fo(t) must be extrapolated 
to t = Aictr = Tijr — m^i which is well above the endpoint 
of the physical region in t for K^ decays. A parameteriza- 
tion that accounts for even higher-order terms is therefore 
desirable. 

The parameterizations belonging to Class I circumvent 
these problems by incorporating additional physical con- 
straints to reduce the number of independent parameters. 
A typical example is the pole parameterization 



nt(t) 



Ks 



K.,s 



f 



(29) 



where the dominance of a single resonance is assumed, and 
the corresponding pole mass My,s is the only free param- 
eter. While for the vector form factor, a pole parameter- 
ization with the dominance of the A"* (892) (M v ~ 892 
McV) is in good agreement with the data, for the scalar 
form factor, there is no such obvious dominance. 
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The most interesting and well-motivated parameteri- 
zations of Class I are those based on dispersion relations. 
These are based on the observation that the vector and 
scalar form factors are analytic functions in the complex 
i-plane, except for a cut along the positive real axis for 
t > tiim — (oik + "t-tt) 2 , where they develop discontinu- 
ities. One can therefore write 



/+,o(*) 



ds (s> 



, Im /+,oQ') 
t-ie) 



subtractions, (30) 



where the imaginary part, Im/ +! o(s'), can be determined 
from data on Kir scattering, and the ultraviolet compo- 
nent of the integral is absorbed into the (polynomial) sub- 
traction terms. In the vector case, the dispersive parame- 
terization turns out to be numerically very similar to the 
pole parameterization due to dominant contribution to 
Im/+(s') from the K*(892). On the other hand, the dis- 
persive parameterization is particularly useful in the scalar 
case, where there is no dominant one-particle intermediate 
state. 

In addition to the analyticity constraints, the scalar 
form factor must satisfy an additional theoretical con- 
straint dictated by chiral symmetry. The Callan-Treiman 
(CT) theorem [73] implies that the scalar form factor at 
t = Ak-k = rn K ~ m w i s determined in terms of Jk/J-k 



and /+(0) up to 0(m Ui d) corrections: 



G = /o(Ak"w) — 



/. 



A' 



U /+(o) 



A. 



CT- 



(31) 



The quantity Act = 0(m u ,d/^f-K) can be evaluated in 
ChPT. At NLO in the isospin limit [44], 



Act = (-3.5 ±8) x 10" 



(32) 



where the error is a conservative estimate of the higher- 
order corrections [73]. Results consistent with Eq. (|32l) 
from NNLO estimates beyond the isospin limit have been 
presented in Ref. 1491741 As discussed in Sect. l3.5.3[ Eq. (13T1) 
provides a useful test of the consistency of the lattice re- 
sults for Jk/ '/tt and /+(0) with experimental data on the 
scalar form factors. 



2.2.5 Dispersive parameterization for the form factors 

Motivated by the existence of the CT theorem, a particu- 
larly appealing dispersive parameterization for the scalar 
form factor has been proposed [71] . Two subtractions are 
performed, one at t = 0, where by definition /o(0) = 1, 
and the other at the CT point, t = Ak-k- Assuming that 
the scalar form factor has no zeroes, this leads to 



f disp/,X _ 



/o (*) = ex P 



t 



A 



Kt, 



-(lnG-G(i)) 



with 



G(t) = 



Ak-k(Aktt — t) 



TV 

3 ds 



<M S ) 



(33) 



(34) 



(s - A Kn )(s -t-ie)' 



With this parameterization, the only free parameter to be 
determined from data is G. 

The phase 4>q(s) can be identified in the elastic region 
with the S-wave (Kn)i = i/2 scattering phase: performing 
two subtractions minimizes the contributions from the un- 
known high-energy phase, which in Ref. 71 is simply and 
conservatively assumed to lie within the interval [0, 27r). 
The resulting function G(t) in Eq. (|3"4"]l does not exceed 
20% of the expected value of In G, while the theoretical 
uncertainties are at most 10% of the value of G(t) |71j . 
The expressions for the leading slope parameters in the 
Taylor expansion as functions of In C are [75] 



A^ = -^[lnC-G(0)] 






(K 



I \2 



2-^G'(0), 

™4 



K' = (KT 6^G'(0)A^ - 3^G"(0) 



IAti 



A 



K- 



where 



IAti 



G(0) = 0.0398(44), 
■G'(0) =4.16(56) x 10" 



Aktt 



G"(0) = 2.72(21) x 10" 



(35a) 
(35b) 
(35c) 

(36a) 
(36b) 

(36c) 



A dispersive representation for the vector form factor 
can be been built in a similar way |75| . Since there is 
no equivalent of the CT theorem in this case, the two 
subtractions are both performed at t = 0. The expression 
analogous to Eq. (|33l) for the vector form factor is 



r sp (i)=exp 



(A + +H(t)) 



with 



H{t) 



m\t [°° ds (f> + {s) 



(s — t — ie) 



(37) 



(38) 



Here the fit parameter is A + = m%df+(t)/dt\t=o and the 
phase 0+(s) is derived from P-wave (Kir)i = i/ 2 elastic 
scattering. As in the case of the scalar form factor, the 
uncertainty on H (t) has a small influence on the determi- 
nation of A+ . The expressions for the leading slopes in the 
Taylor expansion as functions of A+ are |75] . 



A'+ = A+, 
X'i = (A' + 
A"' = (A' + 



2mlH'(0), 
6mlH'(0)X' + 



3m^"(0), 



where 



2m 2 iJ'(0) = 5.79(97) x 10" 4 , 
3m 4 7r H"(0) = 2.99(21) x 10~ 5 . 



(39a) 
(39b) 
(39c) 

(40a) 
(40b) 



The principal results presented in the following sec- 
tions are based on the dispersive parameterizations of 
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Eqs. (|55)) and (|37p . To evaluate the integrals Iki from ex- 
perimental measurements of the parameters InC and A+, 
we use the polynomial expansion given in Appendix [X] A 
detailed discussion of the theoretical uncertainties enter- 
ing into the dispersive parameterization via the functions 
G and H can be found in Refs.lTTl and 1751 

A final remark concerns the isospin-breaking and elec- 
tromagnetic corrections. Throughout this work, a univer- 
sal (isospin-invariant) t dependence is assumed for the 
normalized form factors in the absence of electromagnetic 
effects (i.e., we neglect strong isospin-breaking effects in 
the slopes). Conventionally, the masses appearing in the 
dispersion integrals are chosen to be the charged pion and 
neutral kaon masses, while the correct physical masses are 
used when evaluating the phase space integrals. In princi- 
ple this is not fully correct. For instance, a different cor- 
rection Act should be applied for neutral and charged de- 
cays. However, at present, strong isospin-breaking in the 
slopes can be neglected to well within the experimental 



2.3 Khz and Kt2 beyond the Standard Model 

2.3.1 Effective Lagrangian for semileptonic decays 

The implications of the precision data on Kii and K^ 
decays for SM extensions are most conveniently stud- 
ied within a model-independent effective-theory approach. 
Within this framework, the most general set of weak-scale 
dimension-six local operators contributing to the charged- 
current semileptonic transitions is identified in Ref. 1761 
under the assumptions that the SU(2) x U(l) symme- 
try is linearly realized, and that in the neutrino sector, 
only left-handed neutrinos appear as weak-scale degrees 
of freedom. The resulting low-scale (/x ~ 0(1 GeV)) ef- 
fective Lagrangian governing the semileptonic transitions 
d? — > u 1 £~ i>( for a given lepton flavor £ involves five op- 
erator structures and reads: 



£-di^ui = — 2v2 G F Vij 

(1 + [vi,]ij) ZlIiiVil u l L ^ l d 3 L 

+ [sLJij Irvil u R d 3 L + [an\ij Irvil u L d R 



(41) 



[ti,]ij tRG^viL u R a^ u d 3 L 



h.c. 



where G F /\/2 = g 2 /(8m 2 v ). The effective couplings Vl,r, 
sl,r, and ti, encode information on interactions beyond 
the SM and are of order v 2 /A 2 , where v is the SM Higgs 
vacuum expectation value and A is the new physics scale. 
The coupling vl receives contributions from three gauge- 
invariant weak-scale operators (gauge boson-quark vertex 
correction, gauge boson-lepton vertex correction and con- 
tact four-fermion) while the other couplings are in one- 
to-one correspondence with gauge-invariant four-fermion 
operators at the weak scale. 



In general, the effective couplings carry flavor indices 
and the operators considered here contribute to flavor- 
changing neutral current (FCNC) processes (this is made 
explicit by writing the operators in SU(2)l gauge-invariant 
form at the weak scale). In order to avoid the strong con- 
straints from FCNC, it is convenient to classify the oper- 
ators according to their behavior under the U(3) 5 flavor 
symmetry of the SM gauge LagrangiarQ and organize the 
discussion in terms of perturbations around the U(3) 5 - 
symmetric limit. In practice, we will assume that the un- 
derlying TeV scale new physics has an approximate U(3) 5 
invariance. This can be achieved if flavor breaking is sup- 
pressed by a mechanism such as Minimal Flavor Violation 
(MFV) [77HEI] or by the hierarchy yl flavor > ITeV. 



2.3.2 Phenomenology in the U(3) 5 limit 

We start our discussion by assuming dominance of the 
U(3) 5 -invariant operators. These are not constrained by 
FCNC and can have a relatively low effective scale A. 
Moreover, these operators contribute to a number of pre- 
cision electroweak tests. Therefore, with this analysis one 
can assess the interplay and relative strength of low-energy 
charged-current processes versus other observables (from 
low-energy to the Z pole). 

In the U(3) 5 limit, the phenomenology of charged- 
current processes is greatly simplified [73]: only the SM 
operators (proportional to l + v£) survive in the effective 
Lagrangians both for semileptonic decays (Eq. (J4TJ) ) and 
for muon decay, which process is used to determine the 
Fermi constant. The effects of new physics can therefore 
be encoded into shifts in the values of the effective Fermi 
constants for semileptonic and muon decay: 



Gf = (G F )(°) (1 
G£ = (G F )(°) (1 



vl) 
vl) 



(42) 
(43) 



allowing for different values of vl in semileptonic and 
muon decays. The values of the CKM elements VR cnom as 
determined from semileptonic decays arc affected by both 
these shifts, because semileptonic transitions are normal- 
ized to the Fermi constant G F as determined from muon 
decay. In fact, one has 



T^phcnom t r 






= Vij (1 + v L -vl) 



(44) 



In the U(3) 5 limit, then, a common shift affects all of 
the Vij (as determined from all channels: vector, axial, 
etc.), and the only way to expose contributions from new 
physics is to construct universality tests in which the ab- 
solute normalization of the Vij matters. For light-quark 
transitions, this involves checking that the first row of the 
CKM matrix is a vector of unit length. Therefore, one is 
led to define: 



A 



CKM 



IK 



phenomi2 



\vt 



cnom 1 2 



IK 



phenom 1 2 



1, (45) 



4 I.e., the freedom to perform U(3) transformations in family 
space for each of the five fermionic gauge multiplets: Ql = 



(«i,di), ur, da, Ll — (uljCl), e-R- 
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in terms of the V? cnom determined from semileptonic 
transitions using the standard procedure. The new-physics 
contributions to Z\ckm involve four weak-scale gauge- 
invariant local operators (ip denotes the SM Higgs dou- 
blet), 



Of = ±(L L r<r a L L )(L L ^a a L L ), 

of = (L L ra a L L )(Q L ^a a Q L ), 
of =i{^D»a a V )(L Ll ^ a L L ) +h.c, 
Of = i(^D»a a ip)(Q Llfl a a Q L ) +h.c, 



(46a) 
(46b) 
(46c) 
(46d) 



describing contact four-fermion interactions and gauge 
boson-fermion vertex corrections. Defining 6c\ — r/iV 2 /A 2 
(with rji = ±1), one has 



AcKM = l[af-af-af + af 



Hq 



V 



(47) 



In specific SM extensions, the a.% are functions of the un- 
derlying parameters. Therefore, via the above relation one 
can obtain the constraints from quark-lepton universality 
tests on any weakly-coupled SM extension. 

Each of the a coefficients also contributes to other low- 
and high-energy precision electroweak observables [82] . 
Therefore, we can now address in a model-independent 
way concrete questions such as the following: 

— What is the maximal deviation |Z\ckm| allowed once 
all the precision electroweak constraints have been 
taken into account? 

— Which observables provide the strongest constraints 
on the operators contributing to Z\ckm? 

In order to quantify the significance of the experimental 
constraints on CKM unitarity, we first calculate the range 
of AcKM^cti) allowed by existing bounds from all precision 
electroweak measurements [75]. In terms of the best-fit 
values and the covariance matrix of the &i [82) obtained 
from the fit to electroweak precision data, at 90% C.L. one 
has 

- 9.5 x 1(T 3 < Z\ C km < 0.1 x 10~ 3 . (48) 

This result implies that a deviation from CKM unitarity at 
the level of —1% is not ruled out by precision electroweak 
tests. A closer scrutiny of the precision data shows that 
the blame for large deviations of Z\ckm from zero could 

(3) 

be attributed almost entirely to the operator 0\ , which 
is constrained relatively poorly from LEP2 hadronic cross 
section data, while the other three operators are severely 

(3) • (3) 

constrained: O a by the Fermi constant, V q by hadronic 

(3) 

Z decays, and OV by leptonic Z decays. 

The above discussion implies that even a percent-level 
test of CKM unitarity would provide information not avail- 
able from other precision tests at low and high energies. 
Indeed, by Eq. (|4"7)l. a test of CKM unitarity to better 
than one part in 10 3 (e.g., with a 0.5% determination of 
\V US \ from kaon decays, combined with the 0.02% deter- 
mination of \V u d\ from nuclear beta decays) would probe 
new- physics effective scales A on the order of 10 TeV. As 



we will show in Sect. SJ the current level of theoretical 
and experimental precision in the determination of \V ua \ 
allows this prospect to be realized. 



2.3.3 Beyond U(3) 5 

Corrections to the U(3) 5 limit can be introduced both 
within MFV and via generic flavor structures. In charged- 
current processes, the main effect of U(3) 5 breaking is to 
turn on the chirality flipping (pseudo)scalar and tensor 
structures in Eq. (|4"Tj) . In MFV, the coefficients parame- 
terizing deviations from U(3) 5 are highly suppressed (the 
chirality flip is associated with insertions of Yukawa ma- 
trices). However, such suppression can be compensated 
by a corresponding tan/3 enhancement in models with 
two Higgs doublets of Type II, such as the Higgs sec- 
tor of the minimal supersymmctric extension of the SM 
(MSSM)(3 In this case, the leading non-standard contri- 
bution involves charged Higgs exchange. To one loop, this 
generates the coefficient 



tarr 



[ S R\ 



mpm s 



(1 + eotan/3) m 2 Hi 



(49) 



where eo is a correction factor which is negligible in the 
non-supersymmctric case, while it is 0(l/167r 2 ) in the 
MSSM [S3] (see also [HUES]). In K i2 decay, the H ± - 
exchange amplitude destructively interferes with the SM 
W amplitude; for large values of tan/3 (~ 50) and low 
values of m H ± (~100 GeV), the K& rate can be decreased 
by as much as 5%. 

A highly sensitive probe of U(3) 5 violating structures 
is therefore provided by comparing the value of \V US \ de- 
termined using K^2 decays, which are helicity suppressed, 
and Kgz decays, which are helicity allowedQln practice, to 
minimize the impact of the uncertainties from fx and the 
electromagnetic corrections for K^i it is more convenient 
to consider the ratio 



Rfi23 — 



fl</fir 

/+(0) 



v ud 



[k\ |Krf|o+->o+ 

/J„ 2 [|k.|/+(o)]« 



, (50) 



which makes explicit contact with the results of Sects. 14.11 
and 14.41 The hadronic uncertainties enter through the 
combination (fx/f^/f+iO) and could be reduced if this 
quantity were to be directly computed on the lattice. 
Within the SM i? M 23 = Rnz — 1, while the inclusion of 
Higgs-mediated scalar currents leads to 



R 



fi23 



l - 



*K+ 



tan 2 /3 



1 + eg tan/3 



m. 



h- 



(51) 



5 In such models, tan/3 = V2/V1 is the ratio of the two Higgs 
vacuum expectation values. 

6 In principle, the H exchange amplitude affects K e 3 and 
iiT M 3 decays differently. For K e 3, the effect is totally negligible, 
while for K^, it is substantially smaller than for K^, but not 
totally negligible. However, this effect is well below the present 
theoretical and experimental errors. 
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Note that R^23 also provides interesting constraints on the 
existence of charged right-handed currents that appear at 
NLO in a not-quite-decoupling Higgsless effective theory 

mum. 



3 Experimental inputs 

3.1 Notes on fits to branching-ratio and lifetime 
measurements 

We perform separate fits to world data on the BRs and 
lifetimes for the Kl, K$ and K^. The inputs to our fits 
are the observables actually measured by each single ex- 
periment, such as absolute BRs, ratios of BRs, lifetimes, or 
partial widths. For uncorrelated measurements with sta- 
tistical and systematic errors quoted separately, we add 
the errors in quadrature. In many cases, the results for dif- 
ferent quantities measured by the same experiment have 
correlated errors. The errors are then described by the 
covariance matrix, which must be provided by the ex- 
periment. The free parameters in our fits are the domi- 
nant BRs and the lifetime. In each case, the BRs are con- 
strained to sum to unity. The penalty method [1] is used 
to implement this constraint. Once a first fit has been per- 
formed, scale factors are calculated and used as described 
in the general introduction to the Particle Data Group 
(PDG) compilation [57] ■ The present versions of our fits 
make use of only published measurements. Moreover, for a 
measurement to be included in one of our fits, we require 
the following information to be available: 

— for BR measurements, an explicit discussion of a sat- 
isfactory treatment of radiative corrections, especially 
when they are of comparable size to the experimental 
uncertainty; 

— some details about the estimation of the systematic 
uncertainties; 

— for correlated measurements from the same experi- 
ment, the complete covariance matrix. 



3.2 Dominant Kl branching ratios and tk l 

Numerous measurements of the principal Kl BRs, or of 
various ratios of these BRs, have been published in recent 
years. 

The KTeV experiment has measured five ratios 
of the partial widths for the six main Kl decays 
from independent samples of 10 5 -10 6 events collected 
with a single trigger [91]. The KTeV results for the 
ratios BR(K, l3 )/BR(K e3 ), BR(7r + 7i-7r )/BR(tf e3 ), 
BR(37r°)/BR(ii' e 3), BR(7r+7r-)/BR( J fC e3 ), and 

BR(27T°)/BR(37r ), with total uncertainties ranging 
from 0.4 to 1.1%, are listed in Table [U The six decay 
modes involved account for more than 99.9% of the Kl 
width, so KTeV combines the ratios to determine the 
absolute BRs. We instead use the five measured ratios in 
our global fit to Kl BRs and lifetime. The correlations 
between the errors are provided by the experiment, and 
are taken into account in our fit. 



NA48 has measured the ratio of the K e3 branching 
ratio to that for all Kl decays to final states with two 
tracks [52]. Using a sample of 8 x 10 7 events, they obtain 
BR(AT e3 )/BR(2 tracks) = 0.4978(35). 

Using a sample of 13 x 10 6 4> — >• KlKs decays in which 
the Ks decays to tt + it~ , providing normalization, KLOE 
has directly measured the BRs for the four main Kl de- 
cay channels [3D]. The results depend on the Kl lifetime 
through the geometrical acceptance of the apparatus. The 
values listed in Table[2]were obtained using t k = 51.54 ns 
as a reference value, and depend on the actual value of the 
lifetime as d BR/BR = 0.67 cItk l /tk l ■ KLOE also reports 
results obtained using a constraint on the sum of the four 
measured BRs to solve for tk l , which significantly reduces 
the uncertainties on the BR measurements. In our Kl fit, 
we use the values listed in Table [5] and explicitly include 
the lifetime dependence and other experimental correla- 
tions; we do not make use of the value of tk l obtained 
in [3D]. 

KLOE has also measured tk l directly, by fitting the 
proper decay time distribution using Kl — > 37r° events, 
for which the reconstruction efficiency is high and uniform 
over a fiducial volume of ^0A/3~/ctk l . They obtain tk l — 
50.92(30) ns [88]. We use this value in the fit. 

There are two recent measurements of 
BR(7T+7r-)/BR(i^ 3 ), in addition to the KTeV mea- 
surement of BR(7r + 7r~)/BR(AT e3 ) discussed above. The 
KLOE collaboration obtains BR(7t + 7t-)/BR(A: m3 ) = 
7.275(68) x 10~ 3 [93], while NA48 obtains 
BR(7r+7r-)/BR(A' e3 ) = 4.826(27) x 10~ 3 [SJ. All of these 
measurements are fully inclusive of inner bremsstrahlung. 
The KLOE measurement is fully inclusive of the direct- 
emission (DE) component, DE contributes negligibly to 
the KTeV measurement, and a residual DE contribution 
of 0.19% has been subtracted from the NA48 value to 
obtain the number quoted above. 

We fit the 12 recent measurements listed above, to- 
gether with nine additional ratios of the BRs for subdom- 
inant decays. The complete list of 21 inputs is given in 
Table [2] As free parameters, our fit has the seven largest 
K L BRs (those to K e3 , K&, 3tt°, tt+tt^tt , tt+tt", ttV 
and 77) and the Kl lifetime. Our definition of BR(tt + tt~) 
is now fully inclusive of inner bremsstrahlung (IB), but ex- 
clusive of the DE component. The fit also includes two ad- 
ditional parameters necessary for the treatment of the DE 
component in the radiation-inclusive 7r + 7r _ decay width: 
BR(7r + 7r~7) and BR(7t + 7t - 7de), the branching ratios for 
decays to states with a photon with decay-frame energy 
E* > 20 MeV, and with a photon from DE with E* > 
20 MeV, respectively. Other parameterizations are possi- 
ble, but this one most closely represents the input data set 
and conforms to recent PDG usage. With 21 input mea- 
surements, 10 free parameters, and the constraint that the 
sum of the BRs (except for BR(7r + 7r~7), which is entirely 
included in the sum of BR(7r + 7r~) and BR(7r + 7r _ 7DE)) 
equal unity, we have 12 degrees of freedom. The fit results 
are summarized in Table [3] The fit gives x 2 /ndf = 19.8/12 
(P = 7.1%). The normalized residuals with respect to the 
result of the fit (the pulls) for each input are listed in Ta- 
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Parameter 


Value 


Source 


Ref. 


Pull 


T K L 


50.92(30) ns 


KLOE 


88 


-0.8 


r K L 


51.54(44) ns 


Vosburgh 


89 


+0.9 


BR(K e3 ) 


0.4049(21) 


KLOE 


90 


-1.3 


BR(K^) 


0.2726(16) 


KLOE 


90 


+0.5 


BR( J R- M3 )/BR(iC e 3) 


0.6640(26) 


KTeV 


91. 


-1.1 


BR(3tt°) 


0.2018(24) 


KLOE 


90. 


+2.4 


BR(37T°)/BR(ii'e3) 


0.4782(55) 


KTeV 


91. 


-0.5 


BR(7r+7r-7r°) 


0.1276(15) 


KLOE 


90 


+1.0 


BR(Tv + Tv'Tv )/BR(K e3 ) 


0.3078(18) 


KTeV 


91. 


-0.8 


BR(7r+7r-)/BR(A' e 3) 


0.004856(29) 


KTeV 


91 


+0.3 


BR(7r+7r-)/BR( J R" e 3) 


0.004826(27) 


NA48 


92. 


-0.8 


BR(7r+7r-)/BR(^ M 3) 


0.007275(68) 


KLOE 


93 


-1.5 


BR(ii" e 3)/BR(2 tracks) 


0.4978(35) 


NA48 


94. 


-0.8 


BR(7r°7r°)/BR(37r ) 


0.004446(25) 


KTeV 


91 


+0.6 


BR(7r°7r°)/BR(7r+7r-) 


0.4391(13) 


PDG ETAFIT 


87 


-0.5 


BR( 77 )/BR(3tt ) 


0.00279(3) 


KLOE 


95. 


-0.5 


BR( 77 )/BR(3tt ) 


0.00281(2) 


NA48 


96 


+0.3 


BR(7r + 7r- 7 )/BR(7r+7r-) 


0.0208(3) 


KTeV 


97 


0.0 


BR(7r + 7r- 7 DE)/BR(7r + 7r- 7 ) 


0.689(21) 


KTeV 


98 


+0.2 


BR(7r + 7r" 7 DE)/BR(7r + 7r~ 7 ) 


0.683(11) 


KTeV 


97 


-0.1 


BR(7r + 7r" 7 DE)/BR(7r + 7r- 7 ) 


0.685(41) 


E731 


99 


0.0 



Table 2. Input data used for the fit to Kl BRs and lifetime. For each measurement, the normalized residual with respect to 
the results of the fit is listed in the last column. 



Parameter 


Value 


S 






Corrclal 


ion matrix (%) 






BR(A- e3 ) 


0.4056(9) 


1.3 


-29 


-45 


-30 


+6 


+ 10 


-27 


-27 


+8 


+15 


BR(tf„ 3 ) 


0.2704(10) 


1.5 




-50 





-3 


-10 


-32 


-35 


+13 


-16 


BR(3tt°) 


0.1952(9) 


1.2 






-37 


-1 


+9 


+56 


+63 


-13 


+12 


BR(7r+7r-7r°) 


0.1254(6) 


1.3 








-4 


-16 


-15 


-21 


-5 


-20 


BR(7t + tO 


1.967(7) x 10" 3 


1.1 










+ 14 


+34 


+1 


-3 


+19 


BR(7r+7r~ 7 ) 


4.15(9) x 10~ 5 


1.6 












+16 


+8 


-3 


+74 


BR(7r + 7T- 7 DE) 


2.84(8) x 10" 5 


1.3 














+35 


-10 


+22 


BR(2tt°) 


8.65(4) x 10" 4 


1.4 
















-8 


+10 


BR(ry) 


5.47(4) x 10" 4 


1.1 


















-4 


T K L 


51.16(21) ns 


1.1 





















Table 3. Results of fit to Kl BRs and lifetime. 
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~l I T" 



~i i i I i r~ 



BR(3ti u ) [%] 

~i i I i i i i I r~ 
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Fig. 3. Evolution of average values for main Kl BRs. 



ble [H The evolution of the average values of the BRs for 
K13 decays and for the important normalization channels 
is shown in Fig. [3] 

As seen from Fig. |3l both our fit and the recent PDG 
fit [57] differ substantially in their results from the 2004 
PDG fit [TOO] . This is due to the addition to the world 
data set after 2004 of most of the recent data discussed 
above, and to the subsequent elimination of numerous old 



measurements — many with questionable radiative correc- 
tions and/or unreported correlations — used in previous 
fits. Between 2004 and today, the world-average BRs for 
the K e 3, 3-7T , and ir + Tr~ decays have shifted by +6 a, 
—6 a, and —5 cr, respectively, leading to the following con- 
sequences: 
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— The world-average value for the ratio r^K^/r^^) 
changes from 0.701(8) to 0.6668(28), procuring better 
agreement with the value expected from lepton univer- 
sality, as discussed in Sect. 14.31 

— The world-average value for the amplitude ratio \r)-\ | , 

a parameter of CP violation in the KgK^ system, 
changes from 2.286(18) x 10~ 3 to 2.231(10) x 10" 3 , 
a — 2.7<7 difference. (For this comparison we are inter- 
ested mainly in the effect of the new Kl data, so we 
use the results of the fit to Kg rate data in Table [4] to 
evaluate \rj-\ | for both years.) 

Figure[3]also shows that differences between the results 
of our fit and the most recent PDG fit are minor. These 
fits differ principally in the following respects: 

— The PDG fit uses five of the six KTeV values for the 
dominant Kl BRs, while our fit uses the five ratios 
directly. 

— The PDG fit uses three of the four KLOE values for the 
dominant Kl BRs obtained after imposing the con- 
straint on the sum, as well as the KLOE value of tk l 
obtained by imposing this constraint. We use the orig- 
inal KLOE values with their lifetime dependence and 
omit the value of tk l obtained from the constraint on 
the sum. 

— The PDG fit includes BR(e + e~7) as a free parameter, 
and thus makes use of four measurements not used in 
our fit. 

Our treatment of the contribution from DE in the 7r + 7r~7 
decay is the same as that used by the PDG, but we omit 
the measurement of BR(7r + 7r _ 7)/BR(7r + 7r~) from E731 
[55] . because while the photon energy cutoff in the nu- 
merator is well defined (E* > 20 MeV), the requirement 
used in this measurement for events in the normalization 
channel to contain no calorimeter clusters other than those 
from the tt + tt~ leads to difficulties of interpretation. 



3.3 Dominant Ks branching ratios and tk s 

The KLOE collaboration has measured the ratio 
BR(K S -+ irev)/BR(K s -+ tt+tt-) = 10.19(13) x 10~ 3 
with 1.3% precision |101) . making possible an independent 
determination of |T4s|/+(0) to better than 0.7%. In |102] . 
they combine the above measurement with their measure- 
ment BR(K S -+ ir+Tr-)/BR(K s -+ ir°ir°) = 2.2459(54). 
Using the constraint that the Ks BRs must sum to unity 
and assuming the universality of lepton couplings, they 
determine the BRs for the tt + tt" , 7r°7r°, K e ^, and K^ 
decays. 

Our fit is an extension of the analysis in [102 . We per- 
form a fit to the data on the Ks BRs to 7r + 7r~, 7r°7r°, and 
K e s that uses, in addition to the above two measurements: 

- the measurement from NA48, r(Ks — > ireis)/r(K L — > 
nev) — 0.993(34) [103 , where the denominator is ob- 
tained from the results of our Kl fit; 

— the measurements from NA48, tk s — 89.598(70) ps 
HJBI, and KTeV, t Ks = 89.58(13) ps [IDS], both ob- 
tained without the assumption of CPT symmetry; 



Parameter 


Value 


BR(7T + 7T-) 

BR(7r°7r°) 
BR(K e3 ) 

br(a- m3 ) 

TK S 


0.6920(5) 

0.3069(5) 

7.05(8) x 10" 4 

4.69(6) x 10~ 4 

89.59(6) ps 



Correlation matrix (%) 



-100 



+4 


+3 


+0 


-6 


-6 


+0 




+98 


+1 

+ 1 



Table 4. Results of fit to Ks BRs and lifetime. 



the result BR(K,, 3 )/BR(K e3 ) = 0.6655(15) obtained 
from the assumption of universal lepton couplings, the 
values of A + and In C, the parameters of the dispersive 
representation of the vector and scalar form factors, 
obtained from our fit to form factor data discussed 
in Sect. 13.5.21 and the long-distance electromagnetic 
corrections discussed in Sect. 



The free parameters are the four BRs listed above plus 
tk s ■ With six inputs and one constraint (on the sum of the 
BRs) , the fit has one degree of freedom and gives x 2 /ndf = 
0.015/1 (P = 90%). The scale factor S is not different from 
unity for any of the output values. The results of the fit 
are listed in Table |U 



3.4 Dominant K^ branching ratios and t k ± 

Several recent measurements contribute significant new in- 
formation on the rates for the dominant K decays. In 
addition, we have recently carried out a comprehensive, 
critical survey of the K ± data set, which led to the elim- 
ination of numerous older measurements previously used 
in the fit. The input data used in our current fit to deter- 
mine the dominant K BRs and lifetime are summarized 
in Table El 

The 2003 measurement of BR(K+ 3 ) by E865 [116] 
was the first of the recent-generation measurements of 
semileptonic kaon BRs, and gave a value for \V US \ consis- 
tent with unitarity. The quantity actually measured was 
BR(Ar+)/(BR(7r+7r°) + BR(K+ 3 ) + BR(7r+7r°7r )), where 

one 7T° in the final state was required to undergo Dalitz 
decay. (Throughout the remainder of this section, we use 
7T to denote the charged pion when no confusion results.) 
In 2007, the NA48/2 collaboration published 
measurements of the ratios BR(.K" e 3)/BR(7r7r ) and 
BR(K ^3) /BR(irir°) obtained with simultaneous K + and 
K~ beams [11711121) . For each type of Ki 3 decay (i.e., to e 
and (j,), about 50k K + and 30k K~ decays were collected. 
The results for these ratios in Table [5] are correlated with 
p = +0.19 [122] . The dominant uncertainties are from 



sample statistics. 

ISTRA+ has also performed a measurement of 
BR(K~ 3 ) /BR(tttt°) with 0.6% precision [TJ3]. The result, 
however, has not been officially published and is therefore 
not used in our fit. 

KLOE has measured the absolute BRs for the K e3 and 
K^-j, decays |115j . In — > K + K~ events, K + decays into 
fi + v or 7T + 7r° are used to tag a K~ beam, and vice versa. 
KLOE performs four separate measurements for each K(z 
BR, corresponding to the different combinations of kaon 
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Parameter 


Value 


Source 


Ref. 


Pull 


r K ± 


12.422(40) ns 


Koptev* 


EMI 


+0.9 


r K ± 


12.380(16) ns 


Ott 


107 


-0.3 


T K± 


12.443(38) ns 


Fitch 


108. 


+1.5 


T K ± 


12.347(30) ns 


KLOE 


,109 


-1.2 


BR(X m2 ) 


0.6366(17) 


KLOE 


no. 


+1.1 


BR(7r7r°) 


0.2065(9) 


KLOE 


ED 


+0.5 


BR(7r 7 r°)/BR(A' M2 ) 


0.3329(48) 


PS183 


.112. 


+1.7 


BR(7T7r )/BR( J fi' M 2) 


0.3355(57) 


Weissenberg 


113 


+1.9 


BR(7T7r )/BR(^ 2 ) 


0.3277(65) 


Auerbach 


.114 


+0.5 


BR{K e3 ) 


0.04965(53) 


KLOE 


.115 


-2.1 


BR(A- e3 )/BR(7r7r° + K^ + ttttV) 


0.1962(36) 


BNL-E865 


116 


-0.3 


BR{K ei )/BR{im°) 


0.2470(10) 


NA48/2 


HE2 


+0.6 


BRLsr M3 ) 


0.03233(39) 


KLOE 


.115. 


-3.2 


BR(K„ 3 )/BR(Tnr°) 


0.1637(7) 


NA48/2 


ma 


+1.0 


BR(K„ 3 )/BR(K e3 ) 


0.671(11) 


KEK-E246 


.118 


+0.9 


BR(7T7r 7r ) 


0.01763(26) 


KLOE 


.119 


+0.2 


BR(7T7r 7r°)/BR(7r7r7r) 


0.303(9) 


Bisi 


120 


-0.4 



Table 5. Input data used for the fit to K BRs and lifetime. For each measurement, the normalized residual with respect to 
the results of the fit is listed in the last column. The two 1995 values of t k ± from Koptev et al. are averaged with S = 1.6 
before being included in the fit as a single value. 



charge and tagging decay. As in the case of KLOE's abso- 
lute BR measurements for the Kl, there is some depen- 
dence on the lifetime. The final values for BR(if e 3) and 
BR^X^) in Tableware referred to t k ± = 12.385 ns and 
depend on the actual value of the lifetime as d BR/BR = 
— 0.45 dr K ±/r K ±. The uncertainties are correlated, with 
p = +0.63. 

As seen from Table [5] KLOE has also measured the 
absolute BRs for the important normalization channels 
K+ ->■ tt+tt [HE] and K+ -> p+v [TIP] , using K~ -► 
/i~ v decays to tag. Again, our fit takes into account the 
correlation between these values, as well as their depen- 
dence on the If lifetime. 

One of the primary motives for our critical survey of 
if decay rate data was the poor consistency of the avail- 
able lifetime measurements. The 2007 PDG average value 
for t k ±, 12.385(25) ns, is nominally quite precise. How- 
ever, the error is scaled by 2.1, and the confidence level 
for the average is 0.17% p4] . 

Our survey of the older measurements of t k ± led to the 
elimination of the 1967 result from Lobkowicz et al. |125j . 
because the experiment was much more suited for measur- 
ing the difference between r^-+ and t k - than it was for the 
absolute measurement of either lifetime. The stopped- K + 
measurement described in the 1995 paper from Koptev et 
al. [106] makes use of the surface- muon technique, in which 
kaons are produced and stopped in the same target. The 
two results obtained using different target materials are in 
marginal agreement, with S = 1.6. We have been able to 
identify neither a reason for the discrepancy, nor a reason 
for the exclusion of the measurement, and so we include in 
our fit the average of the two results, with the scale factor 
applied to the error to reflect their disagreement. 

In 2008, KLOE published a new measurement of r^-± 
[TU9"] . The new KLOE result in Table [5] is the average 
of separate measurements for i^ 2 -tagged K + and K~ 



decays using two different techniques for each charge. In 
the first technique, the path length of the tagged kaon 
from production to decay was measured using an inclu- 
sive sample. In the second technique, the tagged kaon was 
required to decay to a final state containing a 7r°, and the 
decay time was measured using the photon clusters in the 
calorimeter. 

Our fit for the six largest K ± branching ratios and 
lifetime uses the 17 measurements in Table [5] The fit has 
seven free parameters and one constraint, giving 11 de- 
grees of freedom. We obtain the results in Table[6] The fit 
gives x 2 /ndf = 25.8/11 [P = 0.69%). The comparatively 
low P-value reflects some tension between the KLOE and 
NA48/2 measurements of the K^ branching ratios. 

For comparison, the most recent PDG fit [57] uses 
32 measurements to determine seven BRs (including that 
for the K e 4 decay to Tr°ir°e u) and the lifetime, and has 
X 2 /ndf = 52/25 (P = 0.13%). The main differences be- 
tween our fit and the current PDG fit, besides the treat- 
ment of the A" ± lifetime data described above, lie in the 
elimination of ten older measurements of various BRs, in- 
cluding the five absolute BR measurements of Chiang et 
al. [126] . In that study, absolute measurements of the six 
largest BRs were performed; the values were improved by 
constraining the sum of the six BRs to unity. Neither the 
original measurements nor the covariance matrix are re- 
ported. Moreover, radiative corrections appear not to have 
been applied when obtaining the BR values. 

Both the significant evolution of the average values 
of the Kis BRs over time and the remaining difference 
between our results and those of the PDG are evident in 
Fig.!] 
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Parameter 


Value 


S 


Correlation matrix 


(%) 




BR(J^ 2 ) 


63.47(18)% 


1.3 


-39 -75 -33 -28 


-36 


+12 


BR(7r7r°) 


20.61(8)% 


1.1 


-26 +61 +38 


-13 


-11 


BR(7T7r7r) 


5.73(16)% 


1.2 


-22 -17 


+36 


-5 


BR(tf e3 ) 


5.078(31)% 


1.3 


+47 


-10 


-13 


BR(K^) 


3.359(32)% 


1.9 




-8 


-4 


BR(7r7r°7r°) 


1.757(24)% 


1.0 






-1 


T K ± 


12.384(15) ns 


1.2 









Table 6. Results of fit to K^ BRs and lifetime. 
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Fig. 4. Evolution of average values for main K^ BRs. 



Experiment 


A' + x 10^ 


A+ x 10 a 


P(A+,A+) 


KLOE 


25.5 + 1.8 


1.4 + 0.8 


-0.95 


KTeV 


21.67+1.99 


2.87 + 0.78 


-0.97 


NA48 


28.0 + 2.4 


0.4 + 0.9 


-0.88 


ISTRA+ 


24.85 ± 1.66 


1.92 + 0.94 


-0.95 



Table 7. Quadratic form factor parameters for Kez decays. 
The values from NA48 and ISTRA+ have been converted for 
use with the notation of Eq. ([28]). The values of p from NA48 
and ISTRA+ were communicated privately |122II127] . 



3.5 Measurements of Kg s form factor parameters 



3.5.1 K e s form factor parameters 



KLOE [Tjg, KTeV [122], NA48 [130], and ISTRA+ [T5T] 

have all performed recent measurements of the quadratic 
parameters A' + and A" of the vector form factor for K e 3 
decays (see Eq. ([25]) ). The data are listed in Table [7] and 
represented graphically in Fig. [5] 

Table [S] gives the results of a fit to the Kl and K~ 
data in the first column, and to the Kl data only in the 
second column. With correlations taken into account, both 
fits give good values of x 2 / n df- The significance of the 
quadratic term is 4.2a from the fit to all data, and 3.5er 



r, 4 

O 

X 

+ 



2 - KTeV 



ISTRA+ 



This fit 




Vx10 J 

Fig. 5. Recent measurements of K e s vector form factor pa- 
rameters. The yellow ellipse shows the result of a fit to all 
data. 



from the fit to Kl data only. Including or excluding the 
K~ slopes has little impact on the values of A' + and A" ; in 



16 



M. Antonelli et al.: Evaluation of \V US \ and Standard Model tests from kaon data 



Kl and K data 
4 measurements 



Kl data only 
3 measurements 



X 2 /ndf = 5.3/6 (51%) x 2 /ndf = 4.7/4 (32%) 



A^ x 10 3 


25.1 ±0.9 


24.9 ±1.1 


\'l x 10 3 


1.6 ±0.4 


1.6 ±0.5 


p(a;,a+") 


-0.94 


-0.95 


I(K° e3 ) 


0.15463(21) 


0.15454(29) 


I{K%) 


0.15900(22) 


0.15890(30) 



Table 8. Average of quadratic fit results for K e3 slopes. 



Experiment 



M v (MeV) 



Ref. 



A+ x 10 



KLOE 870 ± 6 ± 7 [128] 

KTeV 881.03 ±7.11 [129] 

NA48 859 ± 18 [T30] 

ISTRA± 861.6 ± 6.5 ± 5.0 [177] 

Average 871 ± 5 

X 2 /ndf = 3.8/3 



25.74 ± 0.53 
25.17 ±0.48 

26.4 ±1.1 

26.20 ± 0.56 

25.70 ± 0.42 

X 2 / ndf = 2.42/3 



I(K^ 3 ) = 0.15480(18) 
I{K%) = 0.15918(19) 



I{K%) = 0.15478(18) 
I{K%) = 0.15924(19) 



Ref. 



130 

[S3 



Table 9. K e3 form factor parameters for the pole and dispersive parameterizations. Values of A+ from KLOE and NA48 were 
obtained from the pole fit results from those experiments. 



particular, the values of the phase space integrals change 
by just 0.06%. The errors on the phase space integrals are 
significantly smaller when the K~ data are included in 
the average. The results of the fit to all data are plotted 
as the yellow ellipse in Fig. [S] 

All four experiments have fitted their K e3 data us- 
ing the pole parameterization of Eq. (|29l) . and obtain 
the values shown in the left part of Table for the pole 
mass My . The average value is My = 871 ± 5 MeV with 
X 2 /ndf = 3.8/3 (P = 28.9%). The individual values are 
quite compatible with each other, and their average is rea- 
sonably close to the mass of the if* (892). The values of 
the integrals I(K® 3 ) and I(K^ 3 ) as evaluated from the pole 
fit results are just 0.11% higher than the values obtained 
from the quadratic fit results. 

The dispersive parameterization of the vector form fac- 
tor of Eq. (|38p is similar to the pole parameterization, but 
better motivated theoretically. KTeV and ISTRA± have 
used this form to fit their K e3 data, while results for KLOE 
and NA48 can be obtained from the corresponding values 
of My via A + — (m 7T ±/My) 2 . These results are listed 
in the right part of Table [9] together with our average. 
The uncertainties on the values from the individual ex- 
periments do not include the contribution arising from the 
representation of the form factor phase in the dispersive 
parameterization (see discussion in Ref. IB"8"]) . This contri- 
bution is common to all experiments, and is included as 
an additional uncertainty of 0.30 x 10~ 3 on our average 
value of A+ , which is propagated in the evaluation of the 
phase space integrals. 



3.5.2 Kg3 form factor parameters 

KLOE, KTeV, NA48, and ISTRA± have all performed 
measurements of the K^ form factor parameters using the 



quadratic and linear parameterizations for the vector and 
scalar form factor parameters, respectively. The values of 
A^_ and A" from K e3 decays may be averaged with the val- 
ues of A' + , A" , and Ao from K^ decays, with parameter 
correlations taken into account. KTeV and KLOE provide 
K e 3~K^3 averages for the slope parameters; we calculate 
the averages for NA48 and 1STRA±. The K^-K^ aver- 
age values for A+, A" , and Ao from each experiment are 
listed in Table [TO] the la confidence contours for each 
pair of parameters are plotted in Fig. [5] The only signifi- 
cant difference between the results presented in Table \W\ 
and Fig. [5] and the corresponding presentation in Ref. H] 
is that the final ISTRA± systematics are now correctly 
included. As has been previously noted [H I135J . the NA48 
result is incompatible with all of the other measurements. 
Our average of the results from all four experiments is 
plotted in Fig. [6] as the cyan ellipse. This average gives 
X 2 /ndf = 48/9 (P = 3 x 10~ 7 ); the errors on X' + and 
A" are scaled by 1.4 in this case, while the error on Ao 
is scaled by 2.2. If instead of the NA48 K e3 -K^ 3 average 
from Table [Tm only the K e3 measurement from NA48 [130] 
is used, much greater consistency is observed. The result- 
ing average is plotted as the yellow ellipse in Fig. [HI this 
average gives x 2 /ndf = 12.1/8 (P = 14.5%). 

Apart from concerns about the consistency of the data 
set, as described in Sect. 12.2.51 the use of a linear param- 
eterization for the scalar form factor is inherently prob- 
lematic; the use of a Class-I parameterization for the form 
factors such as the dispersive parameterization of Eqs. ([3"3"]> 
and (f3"T|) is greatly preferred. Recently KLOE, KTeV, 
NA48, and ISTRA± have obtained results for the param- 
eters (A + , In C) from fits to their K 1 3 data. These results 
are listed in Table [11] and shown graphically in Fig. [7] 
For each experiment, the results for K e3 and K^ decays 
are averaged. As in Tabled in Table QT] the uncertainties 
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Experiment 


A' + x 10 3 


A'_[ x 10 J 


A x 10^ 


p(K,K) 


P(AV,A ) 


p(a^,a ) 


Refs. 


KTeV 


20.64 ± 1.75 


3.20 ±0.69 


13.72 ±1.31 


-0.97 


±0.34 


-0.44 


1129J 


KLOE 


25.6 ±1.8 


1.5 ±0.8 


15.4 ±2.2 


-0.95 


±0.29 


-0.38 


1132] 


NA48 


24.86 ± 1.88 


1.11 ±0.74 


10.25 ±1.05 


-0.93 


±0.38 


-0.51 


|130lll33| 


ISTRA+ 


24.80 ± 1.56 


1.94 ±0.88 


16.71 ±1.20 


-0.94 


±0.34 


-0.44 


131. 134, 127 


Our avg 


25.02 ±1.12 


1.57 ±0.48 


13.34 ± 1.41 


-0.950 


±0.376 


-0.573 




Our avg, no NA48 A" m3 


25.04 ±0.82 


1.57 ±0.36 


15.90 ± 0.79 


-0.942 


±0.234 


-0.349 





Table 10. Quadratic-linear form factor parameters for K13 decays from KLOE, KTeV, NA48, and ISTRA±. For each ex- 
periment, the results for K^ and K^ decays are averaged. The input K e 3 values from NA48 and all values from ISTRA± 
have been converted for use with the notation of Eq. Q28| l ; correlation coefficients for these measurements were communicated 
privately. Our averages of the results from all four experiments, with and without the NA48 K^ data included, are also listed. 



-KLOE -KTeV - NA48 - ISTRA+ 



our fit 



our fit, no NA48 K 



H3 



? 4 




X.'xIO 3 



? 4 





X n x10 3 



Fig. 6. la confidence contours for measurements of A+, A", and Ao from KLOE, KTeV, NA48, and ISTRA±. For each 
experiment, the results for K e 3 and K^z decays are averaged. Our averages of the results from all four experiments, with and 
without the NA48 K^z data included, are also shown. 



Experiment 


A+ x W 


InC 


P 


Rcf. 


KLOE 


25.70 ±0.57 


0.2038 ± 0.0241 


-0.26 


132 


KTeV 


25.09 ± 0.44 


0.1915 ±0.0116 


-0.27 


68 


NA48 


24.60 ±1.47 


0.1354 ± 0.0133 


-0.24 


133 


ISTRA± 


26.13 ±0.52 


0.2084 ± 0.0134 


-0.46 


ma 


Average 


25.66 ±0.41 


0.2004 ± 0.0091 


-0.33 





Table 11. Dispersive form factor parameters for Ktz decays from KLOE, KTeV, NA48, and ISTRA±. For all experiments 
except NA48, the Kes-K^ average is quoted by the experiment. NA48 does not quote an average; the table lists our average of 
their results, which has S — 2.0 for A+. The values from ISTRA± have been converted for use with the notation of Eqs. (|33[) 
and (|37p : the systematic uncertainties on the ISTRA± measurements are derived from systematic studies using the pole fits. 
Our average of the results from all four experiments, excluding the K^ data from NA48, is also listed. 



on the values from the individual experiments do not in- 
clude the contributions arising from the representation of 
the form factor phase in the dispersive parameterization; 
these contributions are common to all experiments and 
are included as additional uncertainties of 0.30 x 10~ 3 on 
our average value of A+ and 0.0040 on our average value 
of In C. 

The use of the dispersive form factor parameteriza- 
tion clearly illustrates the contrast between the K^ result 
from NA48 and those from the other experiments. The re- 



maining measurements are in agreement. Our average of 
the results from all four experiments, excluding only the 
-K^3 data from NA48, is plotted as the yellow ellipse. This 
average gives x 2 /ndf = 5.6/5 (P — 34.4%). By contrast, 
if the K^3 data from NA48 are included, the average gives 
X 2 /ndf = 25.7/6 (P = 0.026%). On this basis, we exclude 
the NA48 K^ form factor results from the averages used 
to calculate the phase space integrals for the evaluation of 
\V US \ and related tests. 
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O 

c 



0.2 - 



0.15 - 



-KTeV 
-KLOE 
-NA48 
- ISTRA+ 

our fit 




lattice QCD 

0.959(5) 



24 



26 

A x10 3 



Fig. 7. la confidence contours for measurements of A+ and 
InC from KLOE, KTeV, NA48, and ISTRA+. For each ex- 
periment, the results for K e3 and K^ decays are averaged. 
Our average of the results from all four experiments, excluding 
the K^ data from NA48, is also shown. A new, preliminary 
KLOE result |136| is shown as the dashed blue ellipse. This is 
not included in the fit. 



Integral 


A + , A + , Ao 


yl+,lnC* 


Rel. diff. 


I(K 3 ) 


0.15457(20) 


0.15476(18) 


+0.12% 


I(K%) 


0.15894(21) 


0.15922(18) 


+0.18% 


I{K%) 


0.10266(20) 


0.10253(16) 


-0.13% 


KK%) 


0.10564(20) 


0.10559(17) 


-0.05% 


p(K e3 ,K^) 


+0.56 


+0.38 





Table 12. Comparison of phase space integrals evaluated from 
our averages of the results of quadratic- linear (A+, A", Ao) and 
dispersive (A+ , In C) fits. 



Table Q2] lists the values of the phase space integrals as 
computed from the results of our averages of the exper- 
imental form factor data using the quadratic-linear and 
dispersive parameterizations. (In the case of the disper- 
sive parameterization, the uncertainties arising from the 
representation of the form factor phase are included in the 
overall uncertainty for each integral.) For both parameter- 
izations, the correlations between the uncertainties on the 
integrals are described by a matrix of the form 




The table lists the values of the correlation coefficient p 
for each parameterization. 

As seen in the table, when evaluated from the disper- 
sive fit results, the integrals tend to be very slightly dif- 
ferent than they are when evaluated from the quadratic- 
linear fit results. Nowhere is this difference greater than 
0.2%. As expected, the K^ integrals are slightly smaller 
when obtained from the dispersive fit results (Sect. l2~!Q|) . 
Given the advantages of the dispersive parameterization, 




KLOE 0.970(25) 

KTeV 0.982(15) 

NA48 1.039(17) 

ISTRA+ 0.966(16) 



Fig. 8. Values for /+(0) from measurements of InC from 
different experiments, using the Callan-Treiman relation and 
SkIU = 1-193(6). The lattice QCD result /+(0) = 0.959(5) is 
also shown. 



we use the integrals calculated from the dispersive fit re- 
sults for the evaluation of W us \ and related tests. 



3.5.3 The K^j, scalar form factor and tests of chiral 
perturbation theory 

Given a value for Jk/ f-n, the Callan-Treiman relation 
(Eq. (13TJ) ) can be used to obtain a value for /+(0) from a 
measurement of In C, providing a test of consistency be- 
tween scalar form factor measurements and lattice calcu- 
lations. Figure [5] shows the values for /+(0) corresponding 
to the measurements of In C from different experiments, 
using the CT relation and f K /U = 1.193(6) (Sect.|2XTJ). 
For this exercise, the uncertainties on each value of In C 
include the common contribution from the parameteri- 
zation of the form factor phase. The lattice QCD value, 
/ + (0) = 0.959(5) (Eq. ([17])). is also shown. The measure- 
ments of InC from KLOE, KTeV, and ISTRA+ give val- 
ues for /+(0) that are essentially consistent with the lat- 
tice estimate, although they tend to be a little larger. The 
NA48 result, on the other hand, gives /+(0) = 1.039(17), 
which is in contrast with the theoretical expectation of 
Fubini and Furlan that /+(0) < 1 |137j . As noted in 
Sect. 13.5.21 we exclude the NA48 K^ form factor mea- 
surements from the averages used to calculate the phase 
space integrals for the evaluation of \V US | and related tests. 
Our resulting world-average values for A+ and In C are 
listed in Table EI The value InC = 0.2004(91), when 
used with faff* = 1.193(6), gives /+(0) = 0.974(12). 
More generally, the experimental data on InC alone give 
C/W/7r)//+(0) = 1.225(14). This result is completely in- 
dependent of any information from lattice estimates. 

Alternatively, one can perform a fit to the world-average 
value of InC, together with the lattice determinations 
Jk/U = 1-193(6) and /+(0) = 0.959(5), using the con- 
straint given by the CT relation. When performing such 
a fit, we make use of a recent preliminary measurement 
from KLOE [136 of the dispersive K^ form factor pa- 
rameters. This measurement, which is illustrated as the 
dashed ellipse in Fig. [7J is used here in place of the pub- 
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Parameter 


Value 


Corr. matrix (%) 


A+ x 10 J 


25.7 ±0.5 


-27 +10 -10 


InC 


0.208(8) 


-38 +38 


/+(0) 


0.961(5) 


+19 


Jk/U 


1.191(6) 





Table 13. Results of a fit to the experimental average value for 
In C (including a preliminary KLOE measurement) and lattice 
results for /+(0) and fa /fir, using the constraint provided by 
the Callan-Treiman theorem. 



0.213 0.214 

1 ' 1 



0.215 0.216 0.217 

1 ' 1 ' 1 



K L e3 



K L u3 



K s e3 



K*e3 



K*u3 



0.214 0.215 



Fig. 9. Comparison of values for |V^ s |/+(0) for all channels. 
Our average is indicated by the yellow band. 



of common corrections and from correlations in the input 
data set (e.g., from the outputs of the fits to BR and life- 
time measurements). The resulting values of |K ls |/+(0) 
are listed in Table IT4"l and illustrated in Fig. [SI The princi- 
pal experimental result of this review is the average value 



|K s |/+(0) = 0.2163(5), 



(52) 



which has an uncertainty of about of 0.2%. The results 
from the five modes are in good agreement; the fit gives 
X 2 /ndf = 0.77/4 (P = 94%). Table EH gives an approxi- 
mate breakdown of the sources contributing to the total 
uncertainty on the determination of |V^ s |/+(0) from each 
mode. The best single determinations of |V^, s |/+(0) are 
from the K e % and K^ modes, with K e3 giving the slightly 
better determination since there is no contribution from 
uncertainties on the parameters of the scalar form factor. 
The limited precision of the value for |V^ s |/+(0) from K e3 
decays of the K$ decays is entirely determined by the ex- 
perimental uncertainty on the corresponding BR, which is 
dominantly statistical. A better measurement of the BR 
for this decay would allow knowledge of |T4s|/+(0) to be 
significantly improved, since tk s is known very precisely. 
The values of |V^ s |/+(0) from charged kaon decays are 
currently limited in precision both by experimental uncer- 
tainties in the corresponding BRs and by the uncertainty 
in the theoretical estimate of Sg V ^-, . 



4.2 Accuracy of isospin-breaking corrections 



lished KLOE measurement in Table [TT1 Although prelimi- 
nary, we make use of this measurement here as we are only 
interested in demonstrating the power of the consistency 
test offered by the CT relation. The average of the dis- 
persive K^s form factor parameters using the new KLOE 
measurement gives A + = (25.78 + 0.40) x 10~ 3 and InC = 
0.2034(86), with p(A+,lnC) = -0.34 and x 2 /ndf = 7.4/5 
(P = 19.2%). The value of In C is then fitted together with 
the lattice inputs as described above to obtain the results 
in Table [H The fit gives x 2 /ndf = 0.78/1 (P = 38%), 
confirming the agreement between the experimental mea- 
surements of InC and the lattice determinations of /+(0) 
and fx/fir- The values of /+(0) and fic/f-K move only 
slightly, and the uncertainties on these values are slightly 
decreased. 



The average values for |T4s|/+(0) can be computed sepa- 
rately for charged and neutral kaon decays. The ratio of 
these values (calculated without applying SU(2)-breaking 
corrections) then gives an experimental estimate of Sg^fi) 
for comparison with the estimate from chiral perturbation 
theory We obtain 



°SU(2),exp 



0.027 + 0.004, 



(53) 



in good agreement with the theoretical estimate of Ref. 

S SV$) = °- 029 ± °- 004 ( Ec i- ©30 



We observe that the 

uncertainty on the theoretical value of ^g^J) contributes 
significantly to the overall uncertainties on the values of 
|V^ s |/+(0) for charged kaon decay modes (see Table [T4")). 



4 Physics results 

4.1 Determination of \V us \f + (0) 

For each of the five decay modes for which rate measure- 
ments exist, we use Eq. © to evaluate |V^ s |/+(0) from 
the decay rate data in Tables [31 HI and H2 the phase space 
integrals from dispersive fits in Table [T2l the long-distance 
radiative corrections in Table [TJ and the SU(2)-breaking 
corrections of [49] discussed in Sect. 12.2.31 We keep track 
of the correlations between the uncertainties on the val- 
ues of |V^ s |/+(0) from different modes arising from the use 



4.3 Lepton universality 

Comparison of the values of |V^ s |/+(0) computed sepa- 
rately for K e s and K^ modes provides a test of lepton 
universality. Specifically, 



lie 



\Vus\U{Q)]%. 

:iK s i/+(o)]f 3j 



oxp 



Pjf M 3 Ie3 (1 + 2 ^Em) 

r Ke3 j„ 3 (i + 24£) 



(54) 



(Here, 6^m an ^ ^em are * ne electromagnetic corrections 
for kaons of the charge state under consideration.) By 
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Mode 


\Vus\f+(0) 


% err 


BR 


T 


A 


Int 


Correlation 


matrix (%) 


Kl — > 7rev 


0.2163(6) 


0.26 


0.09 


0.20 


0.11 


0.06 


+55 


+10 


+3 


Kl — > TYfj,iy 


0.2166(6) 


0.29 


0.15 


0.18 


0.11 


0.08 




+6 


+4 


Ks —¥ itev 


0.2155(13) 


0.61 


0.60 


0.03 


0.11 


0.06 






+ 1 


K ± -> -Rev 


0.2160(11) 


0.52 


0.31 


0.09 


0.40 


0.06 






+73 


K ± — » TV/M' 


0.2158(14) 


0.63 


0.47 


0.08 


0.39 


0.08 








Average 


0.2163(5) 



















Table 14. Values of |V us |/+(0) as determined from each kaon decay mode, with approximate contributions to relative uncertainty 

°SU(2) 



err) from branching ratios (BR), lifetimes (r), combined effect of 5^ and <5ju(2l (^)i anc ' phase space integrals (Int). 



comparison with Eq. 



' fj,e 



is equal to the ratio gij g\ 



with ge the coupling strength at the W -+ lv vertex. In 
the SM, 77,6 = 1. 

Before the advent of the new BR measurements de- 
scribed in Sects. I3~2l and 13.41 the values of |V us |/ + (0) from 
K e 3 and K^ rates were in substantial disagreement. Us- 
ing the K L and K ± BRs from the 2004 edition of the PDG 
compilation |100) (and assuming current values for the A3 
and 6^), we obtain ry, e = 1.013(12) for K decays and 
1.040(13) for A L decays. 

As noted in Sect. 13.21 the new BR measurements pro- 
cure much better agreement. From the entries in Table [T4l 
we calculate r^ e separately for charged and neutral modes 
(including the value of |14 s |/+(0) from Ag — > i^eu de- 
cays, though this has little impact) and obtain 0.998(9) 
and 1.003(5), respectively. The results are compatible; the 
average value is r Me = 1.002(5). As a statement on the 
lepton-flavor universality hypothesis, we note that the sen- 
sitivity of this test approaches that obtained with -k — > 
Iv decays ((t>) = 1.0042(33) [T3S]) and r -> ivv de- 
cays ((r pe ) = 1.000(4) [I3ni)- Alternatively, if the lepton- 
universality hypothesis is assumed to be true, the equiva- 
lence of the values of |V^ s |/+(0) from A e3 and A m3 demon- 
strates that the calculation of the long-distance correc- 
tions 6 em is accurate to the per-mil level. 



0.228 - 



V ud (0 + ^0 + ) 



V us (Kb) 



0.226 - 



0.224 



1 j- fit with 
1 J unitarity 



N^ 



j 1 I 1 1 i_ 



_i 1 i_ 



0.972 



0.974 



0.976 



v ud 
Fig. 10. Results of fits to \V u d\, \V ua \, and |K S /Fud|- 



discussed in Sect. 14.41 together with the lattice estimate 
f K /f„ = 1.193(6) (Sect. [2"XT1) . Thus we have 



I V as I = 0.2254(13) 
\V US /V ud \ = 0.2312(13) 



[K e3 only], 
[K i2 only]. 



(56) 



4.4 Determination of \V US /V ud \ X f K /U 



As noted in Sect. 12.11 Eq. ([2]) allows the ratio |U us /14d| x 
/k/Itt to be determined from experimental information on 
the radiation-inclusive K12 and 71^2 decay rates. The lim- 
iting uncertainty is that from BR(A m2 ( 7 )), which is 0.28% 
as per Table [fl] Using this, together with the value of t k ± 
from the same fit and r(w ± -> n±v) = 38.408(7) /xs _1 [57] 
we obtain 



\V US /V ud \ x f K /U = 0.2758(5). 



(55) 



Finally, we use the evaluation \V u d\ — 0.97425(22) from 



a recent survey [140] of half-life, decay-energy, and BR 
measurements related to 20 superallowed + — > + nu- 
clear beta decays, which includes a number of new, high- 
precision Penning-trap measurements of decay energies, 
as well as the use of recently improved electroweak ra- 
diative corrections [141] and new isospin-breaking correc- 
tions [142] . in addition to other improvements over past 
surveys by the same authors. Our fit to these inputs gives 



\V ud \ = 0.97425(22), 

\V US \= 0.2253(9) [K a ,K e2 ,0^ 



(U 



(57) 



4.5 Test of CKM unitarity 

We determine \V US \ and \V ud \ from a fit to the results 
obtained above. As starting points, we use the value 
|K s |/+(0) = 0.2163(5) given in Table [H together with 
the lattice QCD estimate / + (0) = 0.959(5) (Eq. ([17])). 
We also use the result \V US /V ud \ x f K /f„ = 0.2758(5) 



with x 2 / ndl = 0.014/1 (P = 91%) and negligible corre- 
lation between |V^d| and \V US \. With the current world- 
average value, \V u b\ = 0.00393(36) [57], the first-row uni- 
tarity sum is then A C km = \V ud \ 2 + |K S | 2 + |Kb| 2 - 1 = 
—0.0001(6); the result is in striking agreement with the 
unitarity hypothesis. (Note that the contribution to the 
sum from \V u b\ is essentially negligible.) As an alternate 
expression of this agreement, we may state a value for 



M. Antonelli et al.: Evaluation of |K, S | and Standard Model tests from kaon data 



21 




Excluded by R^j 

■ 95% CL 

■ 67.28% CL 



1.22 



200 



400 

m H+ (GeV) 



Fig. 11. Regions in the (m H ±,tan/3) plane in two-Higgs- 
doublet models excluded by the present result for 7l M 23 • 



Gckm = G^\V ud \ 2 + \V us \ 2 + \V ub \ 2 . We obtain 

Gckm = 1.16633(35) x 10" 5 GeV" 2 , (58) 

with G M = 1.166371(6) x 1(T 5 GeV" 2 [3J. 

It is also possible to perform the fit with the unitar- 
ity constraint included, increasing by one the number of 
degrees of freedom. The constrained fit gives 

\V US \ = sin 9 C = A = 0.2254(6) [with unitarity] (59) 

and x 2 / ndf = 0.024/2 (P = 99%). This result and that 
obtained above without assuming unitarity are both illus- 
trated in Fig. [TUl 

At this point, using Eq. ([4"T[) and the phenomenological 
value Ackm — —0.0001(6), it is possible to set bounds on 
the effective scale of the four operators that parameterize 
new physics contributions to Z\ckm- We obtain 



A > 11 TeV 



(90% C.L.). 



As noted in Ref. [751 for the operators O), , O , , and 



(3) n (3) 



>(3) 



^l 



Olpq (see Eqs. (|46[)). this constraint is at the same level 



as the constraints from Z-pole measurements. For the 

(3) 

four-fermion operator 0\ , Z\ckm improves upon exist- 
ing bounds from LEP2 by an order of magnitude. 

4.6 Bounds on non-helicity-suppressed amplitudes 



As noted in Sect. 12.3.31 an empirical value for the ratio 
R[i23 (Eq. ([50)) ) can be used to exclude regions of the 
(ni H ± , tan /3) parameter space in models with two Higgs 
doublets, such as the minimal supersymmetry extension 
of the SM (Eq. ([ST]) ). Operatively, we evaluate i? M 23 via 
a fit akin to that described in Sect. 14.51 but with sepa- 
rate parameters accounting for the values of \V US \ from 
Ku, and K^ decays. The fit then has three free param- 
eters: the value of \V US \ from Ki$ decays, the value of 



1.2 - 



1.18 - 




Fig. 12. Results of fit for /+(0) and }k/ fir, given experimen- 
tal data and first-row CKM unitarity: (yellow) no input from 
lattice; (blue) /+(0) = 0.959(5) as input to the fit and no input 
for /k//tt; (red) fx/fn = 1.193(6) as input and no input for 
/+(0). The grey bands illustrate the reference values for /+(0) 
and /k/Ztt- 



\V US /V u d\ from K^2 decays, and the value of \V u d\ from 
+ — > + nuclear beta decays. The input values used for 
\V US \ and IVus/K^I are from Eq. ([56")) and include the rel- 
evant lattice constants. The contribution to non-helicity- 
suppressed K13 decays from charged Higgs exchange is 
negligible, so we include as a constraint in the fit the first- 
row unitarity condition on the value of \V US \ from K^ 
decays: \V ud \ 2 + \V us \ 2 Ki3 + \V ub \ 2 = 1. Expressing the re- 
sults of the fit in terms of | V us \ from K g-j, decays and the 
ratio -7^23, we obtain 



|K S |= 0.2254(8) [i^ 3 ,0 H 
i? M 23 = 0.999(7) \K^\. 



+ , unitarity], 



(61) 



(60) The fit gives x 2 /ndf = 0.0003/1 (P = 99%), with p = 
—0.55 between the parameter uncertainties in the stated 
basis. The regions of the (m H ± , tan/3) parameter space in 
models with two Higgs doublets excluded at the lcr and 
95% CLs by this result for i?^23 are shown as the shaded 
area in Fig.[TT] The bound is obtained setting e = 1/167T 2 
in Eq. (f5i"j). as expected in the MSSM. Note that this result 
excludes the region at low m# ± and large tan f3 favoured 
by B^tv (135) . 



4.7 Determination of Standard Model values for 
/_l_(0) and JkIS-k from experimental data 

Equation @, which in the SM relates the ratio of Km 
and tt£2 decay rates to the ratio |V us /V U( j| x fx/fm can 
be rewritten 



in 



(|K s |/+(0)) 5 



1 



V, 



ad | 



/+(o) 



2 
2 X f2 ' 



I 



(62) 



22 



M. Antonelli et al.: Evaluation of \V US \ and Standard Model tests from kaon data 



where Qti is the ratio oi r{K^2(^)) to /"'(tt^m) with phase 
space differences factored out and radiative corrections 
applied. As per the discussion in Sects. 12.11 and 14.41 for 
K n2M and 7T M 2( 7 ) decays, Q^ a = 0.07604(26). 

Equation ([62]) can be used together with the experi- 
mentally derived values for Qp_i from K^ and 71^2 rates, 
\V us \f+(0) from K i3 decays, and \V ud \ from 0+ -> 0+ nu- 
clear beta decays to determine the values of the hadronic 
constants /+(0) and fx/f-w- A straightforward calculation 
gives 

(f K /U)/U(0) = 1-242(4). (63) 

This value essentially depends only on experimental K^ 
and 7r^2 rate data and radiative corrections. CKM uni- 
tarity is not assumed, although the equality of \V US \ in 
Kfja and K& decays is. If we further assume the valid- 
ity of our reference lattice value /+(0) = 0.959(5), we 
obtain fx/fn — 1.192(8), while if we assume the valid- 
ity of our reference value fx/f* — 1.193(6), we obtain 
/+(0) = 0.960(6). 

To go further, we impose first-row CKM unitarity as an 
additional constraint. We can then perform a fit with <3 M 2, 
\V us \f+(0), \V u d\, /+(0), and fx/fn as parameters, and 



two constraints: \V u d\ 



\V US 



|Kf,| 2 = landEq. (EH). 



Up to two of the parameters can be left free for a solution 
with no degrees of freedom. Taking these to be /+(0) and 
fx/U, we obtain 



/+(0)= 0.959(5), 
fx/U = 1-192(6) 



[with unitarity], 



(64) 



with a correlation coefficient of p = +0.84. This solution 
is illustrated as the yellow ellipse in Fig. [121 Once again, 
we note that this result is obtained without any recourse 
to prior estimates of /+(0) or fx/fn- The lattice values of 
these constants are illustrated as the grey bands in Fig. [12] 
for the purposes of comparison only. 

With either of the reference values of /+(0) or fx/fir 
as an additional input, the fit has one degree of freedom. 
With /+(0) = 0.959(5) as input, the fit gives \V ud \ = 
0.97425(18), |K s |/+(0) = 0.2163(4), /+(0) = 0.9594(34), 
and fx/fn = 1-192(5), with p = +0.76 the correlation co- 
efficient between /+(0) and fx/f-x, and x 2 = 3x 10 _4 /1 
(P = 99%). This is illustrated as the blue ellipse in Fig. [12 
With fx/f-Tt = 1-193(6) as input instead, the fit gives 
\V ud \ = 0.97427(17), \V us \f+(Q) = 0.2163(5), / + (0) = 
0.960(4), and f K /U = 1-192(4), with p = +0.75 and 
X 2 = 0.023/1 (P = 88%). This is illustrated as the red 
ellipse in Fig. 1121 In summary, our chosen reference val- 
ues for /+(0) and fx/f-w are a near-perfect match with 
experimental data and SM assumptions. 



5 Conclusions 

At present, the experimental precision of data on leptonic 
and semileptonic kaon decays is nicely matched below the 
percent level by the theoretical precision, allowing us to 
perform very precise measurements of SM parameters and 
to set stringent bounds on beyond-SM physics. In this 



work, we have presented an updated analysis of the fol- 
lowing three main subjects: (i) the overall determination 
of \V US \, with and without imposing CKM unitarity; (ii) 
the comparison between the values of \V US \ as determined 
from K?3 and K ^ decay data, and the corresponding con- 
straints on deviations from the V — A structure; (iii) tests 
of lepton universality in K13 decays. The final results of 
our analysis are presented in detail in Sect. 2] and will 
not be repeated here. As a general conclusion, we sim- 
ply emphasize that the 0(10 TeV) bound on the scale of 
new physics (Eq. (|60p ). which follows from the verification 
of the first-row CKM unitarity condition, represents one 
of the most stringent constraints on physics beyond the 
Standard Model. 

In addition to the three main topics mentioned above, 
our analysis has allowed some interesting cross-checks of 
both analytic and lattice-based theoretical results rele- 
vant to kaon physics, most notably the determination of 
the mass-induced isospin-breaking effect between charged 
and neutral semileptonic kaon decays (Eq. (153]) ) and the 
indirect determination of the kaon form factors shown 
in Fig. [T^l We emphasize that both of these results, 
which agree well with the current theoretical estimates, 
are strictly valid only under a pure SM hypothesis. 
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A Polynomial approximations for the phase 
space integrals Iki 

In this appendix, we give the explicit forms used to eval- 
uate the phase space integrals Ixe using the form factor 
parameters obtained from fits with different parameteri- 
zations. In all cases, we insert the expressions for / +i o(£) 
into the integrand of Eq. (|27[) and obtain a polynomial 
expression in the form factor parameters for the integral 
Ixe, with coefficients depending on the mass values mx 
and mg. 
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The form of the phase space integral used with the 
polynomial parameterization up to cubic terms in t is 



Iki = z + zi\' + + z 2 (X' + + A") + z 3 X' + X'l + z 4 X 
+ z 5 X+ + z 6 X' + X^ + Z7 X'lX'^ + z 8 X'f 
+ zgX a + zio(A + A ) + znA A + Z12A0 
+ 2 i3A + zi4A A +zi5A A + zi6A . 



//2 
+ 



(65) 



The numerical values of the coefficients z are listed in 
Table US] For K e3 decays, we take z 9 . . . zn to be zero. 
We also use the above form with to evaluate phase space 
integrals from the results of fits using the pole parameter- 
ization. In this case, we first perform a Taylor expansion 
of Eq. (|^|) up to cubic terms. 

The form of the phase space integral used with the 
dispersive parameterizations of Eqs. (J33I) and (l37l) is 



Iki = c + c x A+ + c 2 A% + c 3 A 3 , 



c±A% 



+ J + ■ --t (66) 

+ c 5 InC + c 6 (lnC) 2 + c 7 (lnC) 3 + c 8 (lnC) 4 . 

The numerical values of the coefficients c are listed in Ta- 
ble [THJ Again, for K e3 decays, we take C5 . . . eg to be zero. 
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Coeff. 


I{K%) 


I(K l > 3 ) 


I(Ka) 


UK*) 


Coeff. 


HKs) 


HKb) 


zo 


0.14478436 


0.14081706 


0.09345981 


0.09083246 


Z9 


0.21195460 


0.20545016 


Z\ 


0.50059594 


0.48654454 


0.30229496 


0.29365924 


ZlO 


0.37262925 


0.36033462 


Zl 


0.69266614 


0.67217878 


0.48877555 


0.47401926 


Zll 


1.52866642 


1.47403632 


Zz 


2.39194070 


2.31634083 


1.80851699 


1.75020047 


Zl2 


1.73806537 


1.67074681 


Zi 


2.35425967 


2.27427126 


1.84604117 


1.78213612 


Z13 


0.50955547 


0.49134544 


Zh 


0.79731357 


0.77211361 


0.60283900 


0.58340016 


Zl4 


2.31742049 


2.22766241 


Z6 


3.13901290 


3.03236168 


2.46138823 


2.37618149 


Z15 


5.65071063 


5.41411258 


Z7 


6.74124302 


6.49464837 


5.40291925 


5.20184325 


ZlG 


4.82367980 


4.60616108 


z% 


5.13492355 


4.93286555 


4.17536001 


4.00843393 









Table 15. Coefficients for the evaluation of phase space integrals from polynomial form factor parameters, for use with Eq. ([65 



Coeff. 


I(K%) 


I(K u e3 ) 


m%) 


KK*) 


Coeff. 


m%) 


HKs) 


CO 


0.14530227 


0.14125098 


0.09323787 


0.09059706 


C5 


0.01797945 


0.01724105 


Cl 


0.50409859 


0.48947133 


0.30495081 


0.29587728 


C6 


0.00545007 


0.00515380 


C2 


1.39922105 


1.35593685 


0.98846227 


0.95712982 


Cl 


0.00128587 


0.00119853 


C3 


3.22924652 


3.12170948 


2.44345185 


2.36027254 


C8 


0.00025237 


0.00023179 


c 4 


6.36982500 


6.14083977 


6.53169253 


4.81427789 









Table 16. Coefficients for the evaluation of phase space integrals from dispersive form factor parameters, for use with Eq. (|66|) . 
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